Distribution of velocities and acceleration for a particle in Brownian correlated 

disorder: inertial case 
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We study the motion of an elastic object driven in a disordered environment in presence of both 
dissipation and inertia. We consider random forces with the statistics of random walks and reduce 
the problem to a single degree of freedom. It is the extension of the mean field ABBM model 
in presence of an inertial mass m. While the ABBM model can be solved exactly, its extension 
to inertia exhibits complicated history dependence due to oscillations and backward motion. The 
characteristic scales for avalanche motion are studied from numerics and qualitative arguments. To 
make analytical progress we consider two variants which coincide with the original model whenever 
the particle moves only forward. Using a combination of analytical and numerical methods together 
with simulations, we characterize the distributions of instantaneous acceleration and velocity, and 
compare them in these three models. We show that for large driving velocity, all three models share 
the same large-deviation function for positive velocities, which is obtained analytically for small and 
large m, as well as for m — 6/25. The effect of small additional thermal and quantum fluctuations 
can be treated within an approximate method. 



I. INTRODUCTION 

The dynamics of a large class of classical and quan- 
tum systems can be modeled within the description of 
an elastic manifold driven by an applied external force 
through a disordered medium [ll-Q- Some examples are 
domain walls in magnetic systems in the presence of time- 
dependent magnetic fields [J], flux- line lattices in type- 
II superconductors driven by an applied transport cur- 
rent [5| , charge-density waves in solids in an electric field 
0] , pinned or driven Wigner crystals |8l-lll| , disloca- 
tions in metals [13, interface between two fluids in a 
porous medium |13|, e arthquakes [l|, and crack fronts 
in brittle materials [IJ]. In all these systems, the com- 
petition between elastic forces, quenched disorder and 
external driving shapes the dynamics. As a result, the 
response is usually complicated. 

If the driving force is sufficiently small, the system is 
trapped due to disorder in a metastable state. When in- 
creasing the external driving, some weakly pinned parts 
will start moving. They will be stopped by elastic forces 
that describe interactions between weakly and strongly 
pinned regions of the manifold. Further increase of the 
driving usually results in jumps of a segment of the sys- 
tem, and avalanche motion. 

One example occurs in soft magnets. When smoothly 
increasing the magnetic field (H), the magnetization (M) 
changes in an irregular way. This process can be ex- 
plained by considering the motion of domain walls sepa- 
rating regions of opposite magnetization. The derivative 
of the magnetization with respect to the magnetic field, 
dM/dH, is known as Barkhausen noise and can be re- 
lated to avalanche motion [1^, [Ig . An important step to- 
wards modeling the dynamics in these systems was made 
by Alessandro, Beatrice, Bertotti and Montorsi (ABBM). 
On a phenomenological basis they introduced [17|i 1131 ^ 
Langevin equation for the velocity of a single degree of 
freedom, i.e. a particle, which represents the center of 



mass of the domain wall. It is simple enough to allow for 
an exact solution. Their approach is known as ABBM 
model. It was successful in explaining the distribution of 
sizes and the duration of pulses in the Barkhausen signal, 
both for an extremely small and for a finite increase-rate 
of the external field [i. IT9l-[2H . In the ABBM model, the 
probability of the instantaneous domain-wall velocity is 
found to be P(u) ^ u~" cxjp{^u/uo), where a = 1 — f . 
Here v is proportional to the rate of increase of the field 
and Mo is some characteristic cutoff. The avalanche sizes 
S are distributed according to P{S) ~ S'^'^ up to some 
large-scale cutoff, with r = (3 — v)/2. For vanishing 
rate w — > 0"*" , different samples of different materials were 
found to be characterized by universal exponents a and 
T, regardless of the specific microscopic details about the 
sample structure. 

In their phenomenological theory, ABBM assumed 
that the random-force landscape seen by the particle has 
the long-range correlations of a Brownian motion, while 
the original disorder seen by the domain wall is of short- 
ranged nature. This approximation, made in order to 
explain experiments, turns out to be justified in some 
cases. First, in the limit of an interface with infinite- 
ranged interactions (i.e. a fully connected lattice model) 
it was shown that the ABBM model becomes exact, i.e. 
it describes exactly the center-of-mass motion [201 ■ It 
is believed to provide a mean-field model, which should 
be valid in particular to describe domain walls in situa- 
tions where the long-ranged dipolar forces generate long- 
ranged elasticity and puts the system at its upper critical 
dimension [22]. 

Recently, two of us have developed a field theoretic 
approach to describe and compute avalanche-size distri- 
butions and velocity distributions for elastic interfaces of 
internal dimension d in short-ranged disorder [23l - (25| . It 
was shown that in the quasi-static limit w — >■ 0^ the veloc- 
ity of the center of mass is indeed described by the ABBM 
model at and above the critical dimension dc, with cor- 



rections subdominant in the spring constant (parameter 
/i below). Deviations become important for d < dc, and 
were computed in a d ~ dc — e expansion, where dc — 4: 
for short-range elasticity, and dc = 2 in case of dipolar 
forces. The theory also allows to predict the spatial de- 
pendence of avalanches [2J,[2al which cannot be obtained 
from the ABBM model. It also provides an independent 
exact solution of the ABBM model at any driving velocity 
[26| based on the Martin-Siggia-Rose (MSR) formalism 
via the solution to a non-linear saddle point equation, 
called the instanton equation. 

The ABBM model, and the subsequent field theoretical 
approach, provides a good descriptions of avalanche mo- 
tion in classical systems evolving with the simplest over- 
damped dynamics. One would like to extend these the- 
ories to describe elastic systems with a more general dy- 
namics, including inertial and retardation effects, and to 
describe avalanche dynamics in quantum systems. Stud- 
ies of classical models with stress overshoots [23, [2^ has 
shown that although the depinning transition may be 
not too much affected in the thermodynamic limit, the 
avalanche-size distributions can be quite different. 

Retardation effects arc important for instance in mag- 
nets. Apart from universal power laws discussed above 
(characterized by the exponents a and t), pulses of dif- 
ferent durations in Barkhausen noise are expected to col- 
lapse on the same curve after proper rescaling |29| . How- 
ever, in some experiments on ferromagnetic alloys the 
pulse shape is found to be asymmetric J29l - l3ll |. This 
asymmetry was explained to be a transient effect of eddy 
currents [32, ^M ■ Namely, the domain wall motion gener- 
ates eddy currents. The response is not immediate, but 
instead finite-time delays exist after the corresponding 
wall displacement is made. These effects of retardation 
can be taken into account by introduction of a negative 
mass of the domain wall [S^, ISJI • 

Although the above-mentioned effects are important 
in some samples, inertial effects are always present in the 
domain wall dynamics. A domain wall is characterized by 
the so-called Doring mass, which is due to gyromagnetic 
effects [SJ]. However, inertial effects are often neglected 
with respect to a larger damping present in the system, 
and simplified models excluding the mass are studied. In 
many other systems the dynamics is only weakly dissipa- 
tive and inertial effects can be important. Some examples 
are geological faults, motion of contact lines of a droplet 
on a dirty rough surface and crack fronts in brittle ma- 
terials. Domain walls with an internal degree of freedom 
also exhibit a non-trivial dynamics reminiscent of inertial 
effects m. 

The description of avalanches in quantum systems with 
quenched disorder is also a challenge. There is great cur- 
rent interest in non-equilibrium quantum systems and 
their full counting statistics J36l - l40l |. Higher moments 
of the noise have been measured in avalanche processes 
and exhibit some resemblance to their classical counter- 
parts [4ll |. Out of equilibrium elastic quantum systems 
in presence of disorder and a bath have been studied in 



the thermal, and quantum-creep regimes [42, |43| where 
the driving force is small and the dynamics is slow and 
governed by the time scales set by by thermal or quan- 
tum tunneling over barriers. For a fixed driving force 
above the depinning threshold however, there are no bar- 
riers. To study avalanches, it is convenient to drive the 
system with an external spring at fixed but small ve- 
locity V [2j|. Then the effective driving force changes 



in time, and the spring provides a restoring force that 
keeps the system near the depinning transition. In the 
stationary state the system is temporarily pinned, then 
unpins and jumps to the next metastable state, and so 
on. In that situation, while thermal or quantum fluctu- 
ations may help trigger an avalanche (see e.g. p. 319 in 
Refs. [j| or [ij]) they should be less important during 
the avalanche process itself, which usually involves much 
faster time scales than that of barrier crossing. During 
the avalanche the system is rolling down the potential 
hill, with possible overshoots and oscillations due to in- 
ertia. Thus a semi-classical equation of motion keeping 
only inertia and damping into account should be a rea- 
sonable starting point. 

Given these motivations, in this paper we study the 
ABBM model in presence of inertia. We consider the 
motion of a particle representing the center-of-mass po- 
sition of an interface, that is driven by a spring at velocity 
u > in a Brownian-correlated random-force landscape. 
The feature which makes the ABBM model solvable is 
that the motion is always forward. In presence of inertia 
this property is lost. The disorder thus generates non- 
trivial correlations in time when the particle visits the 
same positions several times. To make progress we thus 
consider two variants of the model. 

One variant is a model "on a tree", i.e. such that 
when the particle changes the direction of motion, it ex- 
periences a different Brownian disorder. Although it may 
seem artificial, it could in fact be of relevance for inter- 
faces since different parts of an interface are exposed to 
different disorder potentials, and in presence of inertia 
the backward motion of the center-of-mass does not have 
to involve the same segments of the system as the forward 
one. The advantage of this model is that it can again be 
studied using a Fokker-Planck equation, which however 
does not appear to be exactly solvable. We determine 
the joint distribution of velocity u and acceleration a: (i) 
in perturbation theory at small and large mass, (ii) for 
large driving velocity v, and (iii) numerically. We then 
compare with a numerical solution of the original model, 
i.e. the ABBM model with inertia. 

The second variant we call the yu model. It is the 
model for which the method developed in |24| naturally 
extends. The non-linear instanton equation is now a dif- 
ferential equation of second order in the time variable. It 
is the saddle point equation of the MSR action for the 
ABBM model with inertia under the assumption that 
the particle moves in the direction of the drive only. We 
arc unable to solve it exactly for generic values of the 
mass. We solve it: (i) in perturbation at small and large 



mass, and (ii) for a magic value of the mass where ex- 
act solutions exist, related to the Abel equation. It is 
also easy to solve numerically and from it we obtain, for 
that model, the Fourier-Laplace transform of the veloc- 
ity distribution. Also, we calculate exactly the moments 
characterizing the distribution function (a'^'u") for arbi- 
trary V and mass. The only unpleasant feature of this 
model is that due to backward motion, complex veloci- 
ties appear. As long as they have a small probability, e.g. 
for large v or small mass, it gives the correct physics. In 
fact, by comparing with numerics, we find that this model 
provides quite interesting approximations to the ABBM 
model even for not so small values of mass and velocities. 

Although the three models, namely the original one, 
the tree model and the y/ii model, do correspond to differ- 
ent ways to treat the negative velocities, it appears that 
they share the same large-deviation function at positive 
instantaneous velocity. The latter describes the large v 
limit (i.e. driving velocity) of the probability distribution 
of the (instantaneous) velocity and acceleration. Hence 
we conjecture that we have obtained in this paper the 
exact large-deviation function for the ABBM model with 
inertia in the positive-velocity domain. This conjecture 
is explained and argued for in details, and supported by 
numerics. We find that the large-deviation function is 
determined by the nonlinear instanton equation and we 
obtain its analytical form: (i) in perturbation for small 
m and for large m; (ii) for the magic value of the mass. 
In addition we discuss for all three models the probabil- 
ity that the particle, starting with given acceleration and 
velocity, reaches zero velocity before or at time t . We 
refer to the latter as the exit probability. 

Although our calculation is performed at zero temper- 
ature and h, the distribution of velocities and accelera- 
tions in an avalanche is expected to be robust, and should 
survive at low temperatures, as well as in the presence 
of c^uantum fiuctuations. More precisely, from the above 
discussion, it should be valid as as long avalanche du- 
rations remain small compared to barrier-crossing time- 
scales, which is the regime studied in this paper. A more 
complete theory however, yet to be worked out, would 
need to incorporate several additional effects: (i) the 
renormalization of disorder by fluctuations; (ii) in the 
under-damped limit the total avalanche duration may be 
notably increased as the system oscillates before settling 
into the next metastable state; (iii) the scale dependence 
of these effects. Although the present approach is only a 
first step, it is expected to capture some of the effects of 
inertia in classical and quantum avalanches. In particular 
at the end of the paper we show how to incorporate some 
of the thermal and quantum effects in the moving system 
by studying the y/u model in presence of an additional 
thermal or quantum noise. 

The paper is organized as follows. First, in Sec. |lll 
we rederive the distribution of velocities for the ABBM 
model and make the connection with the MSR formalism 
employed in Ref. |2J]. Then we start analyzing inertial 
effects. In Sec. Illll we introduce the ABBM model with 



inertia and analyze the results of a numerical simulation. 
Then, in Sec. IIVI we consider a particle on "the tree" . 
The model is introduced in Sec. IIV Al In Sec. IIVBI we 
solve the corresponding Fokker-Planck equation pertur- 
batively in the inertia, and determine the corresponding 
probability distribution. The large-w limit is discussed 
in Sec. WC\ In Sec. HVDl we solve the Fokker-Planck 
equation perturbatively in 1/m and give a solution for 
the -\Ai-model in the same limit. In Sec. IIV El we solve 
the Fokker-Planck equation numerically, and finally com- 
pare analytical and numerical results in Sec. IIVFI Then, 
in Sec. IV Al we consider the V^-modcl. We start with 
the definition and its basic properties in Sec. IV Bl Then 
we connect the instanton and Fokker-Planck approaches 
in Sec. IV CI In Sec. IV Dl we calculate exactly the mo- 
ments characterizing the distribution function. We solve 
the instanton equation perturbatively in the mass and 
from that we find a perturbative expansion of the distri- 
bution function in Sec. IV El while in Sec. IV Gl we solve 
it exactly for the "magic" value of the mass. In Sec. IV Fl 
we analyze in more detail one fixed value of the mass. 
In Sec. I VII we introduce and discuss the large-deviation 
function. Its perturbative expansion in small and large 
TO is given in Sec. I VI CI and in Sec. IVIDI respectively, 
while in Sec. I VI El we give its exact result for the magic 
value of the mass. Supplementary material is relegated 
to appendices [Alto III 



II. ABBM MODEL 

Before we start considering inertial and dissipative ef- 
fects together, we first review the ABBM model [13, [3 
that neglects inertia, rederive its velocity distribution, 
and recall the connection to the saddle point (instanton 
equation) approach of Ref. |24|-[26l| . 

We study an elastic interface at zero temperature, 
whose center-of-mass position is given by the equation 
of motion 



i]u{t) = F{u{t)) + ^j? [vt - u{t)] , 



(1) 



where ii = du/dt and F is the disorder force. It is Gaus- 
sian distributed with correlations 



[F{u) - F{u')] 



2a\u 



(2) 



In Eq. ([T]), 1] measures dissipation and v is the driving 
velocity. For the specific realization in magnetic sam- 
ples, the interface describes a domain wall and the term 
^ fi^u models the demagnetizing field generated by free 
magnetic charges on the boundary of the sample [J| . In 
general, this term is a restoring force and /i^ the spring 
constant by which the particle (representing the ccntcr- 
of-mass) is driven. 

For V > 0, the Middleton theorem [45| states that the 
particle always moves forward in the steady state (and 
for all t > if its initial velocity at f = is positive). 



The above equations can be solved via the Fokker-Planck 
equation 

dP{u,t) d 



j{u,t). 



dt dii 

Here j{u,t) is the probabihty current 



(3) 



j{u,t) 



u?v U^u\ „,. X d f ail ^, . . 
t^-^)p{u,t)- — (—P{u,t) 



(4) 



The first term is the contribution from the drift and the 
second one from the diffusion. 

Before proceeding further let us recall the main scales 
for the ABBM model, and introduce the appropriate di- 
mensionless units to be used in this paper. Times will be 
measured in units of the relaxation time of the quadratic 
well, 



/i2 



(5) 



Displacements (i.e. the u direction) will be measured in 
units of: 



S",, = o-/^'' 



(6) 



where 5^^ gives an estimate of the large-size cutoff for the 
distribution of avalanche sizes, as defined in [2J]. Veloci- 
ties will be thus measured in terms of a velocity scale set 
by the disorder. 

With these units of time and space the ABBM model 
contains only one parameter, the driving velocity v :— 
v/v^ in dimensionlcss units. Below, we will mostly use 
these units, keeping the freedom to restore dimensionfull 
units when needed. 

Let us now discuss the steady state of the ABBM 
model. In that case dtP = and j{u) = const. How- 
ever, from the condition that the particle always moves 
forward follows that j{u) = 0. Solving Eq. ^ with this 
constraint, one obtains (in dimensionlcss units): 



P{u) 



T{v) 



-e{u). 



(8) 



It is important to note the dramatically different behav- 
ior of lini^,_j.o P(u) for v < 1 and v > 1. In the former 
it is divergent while in the latter it tends to zero. The 
value z; = 1 separates the regime of intermittent motion, 
where the particle is most of the time at rest, from the 
regime where it moves continuously. 

Now we briefly discuss an alternative way of solving the 
Fokker-Planck equation and make connection with the 
approach introduced in Ref. [2^ based on the instanton 
equation, and further studied in [26| . After performing 
the Laplace transform we find 



dP 



9 



M' 



^- + A^P(A) - ^XdxP + -^ydxP + BT = 0. 



dt 



V 



V 



ij^ 



(9) 



Here P{X) = L P{u)e^"du and the boundary terms 



(BT) arc 
BT 



[j{u,t)e 



Xul 



if' 



[uP{u)e^^] 



(10) 



We ignore for the moment the boundary terms and look 
for a solution of the form 



P{\t) 



- pVZ{\,t) 



(11) 



with Z{0,t) = since P(A = 0,t) = 1. After introducing 
dimensionless quantities A' = Xa/r]^i^, and Z' = Za/rj^^ 
and omitting primes, we find |24| 



dZ{\,t) dZ{\,t) 



dt 



dX 



(A - A^) = A. 



(12) 



This equation admits a time-independent solution 
Z{X,t) = Z(A) withZ(O) = 0: 



Z(A) = -ln(l-A) 



(13) 



Hence we recover the result of [2J] for the steady state. 
Doing the inverse Laplace transform of P(A) = (1 — A)"" 
we obtain Eq. (|5]). Now, we can check that the boundary 
terms indeed vanish for A < 1, i.e. in the domain in which 
Z{X) is defined. 

In addition one can make the connection to the instan- 
ton equation as follows. The equation ([T^ can be solved 
by the method of characteristics: Define a function A(i) 
which obeys the differential equation 



dXjt) 
dt 



Xit)-X'it). 



(14) 



Further define Z{t) := Z{X(t),t). Then the total deriva- 
tive, using ([T2|) is 



dZjt) 
dt 



m- 



(15) 



The equation ([T^ is nothing but the instanton equation 
of Ref. [2J|, with here A(t) = u{t) there. It admits the 
solution 



m = 



X 



A + (1-A)e-* 
with boundary condition A(— oo) = 0. In addition. 



Z{t) 



A(i') dt'. 



(16) 



(17) 



where we have defined A(0) = A. Hence if we express 
Z{X) :— Z{t = 0) as a function of A = A(0) we obtain 
precisely ^. In Ref. I24,i25|| Eqs. ^ and ^ were ob- 
tained not as the Laplace transform of the Fokker-Planck 
equation, but by a completely different route using the 
Martin-Siggia-Rose (MSR) dynamical action. This will 
be explained in more details below in the case where in- 
ertial terms are allowed. 



Eq. ((T2)) is solved for any initial condition Z{\,t 
0) = Zo(A) as mi 



with 



Z{\i) = -ln(l-A+Ae-*) + Zo 
Hence, 
P(A,i) = (l-A + Ae-*)-"Po 



A 



A + (1-A)e* 



A 



A+(1-A)e* 



(18) 



(19) 



Using u"(t) ~ 95i'P(A,i)|A=Oj we obtain the decay to the 
steady state after a change in the driving velocity, v = 
vo + e{t){v - vq): 



u{t)=v{l-e-^) + e-*u{Q) (20) 

'iW' = «(1 - e-'f + 2^e-*(l - e-*) + e-2*^i(0)2", 

where the symbol c denotes connected correlation func- 
tions. This is in agreement with the results of Ref. |2a |. 
where it was obtained within the MSR approach. Note 
that for any t > Q Eq. (HH) behaves as P{\,t) ^ 
A{t){-\Y with A{t) = (1 - e-*)-"Po(-^), hence 
P{u,t) ~ A{t)u'"~^ /T{v) and the current at the origin 
vanishes, which justifies ignoring the boundary terms 
above Hg. 



III. ABBM MODEL WITH INERTIA 

A. Definition of the model 

In this section, we consider the generalization of the 
ABBM model to include the effect of inertia. The equa- 
tions of motion in the laboratory frame are 



du{t) 

dt 
da{t) 



ait), (21) 

^i^[v - u{t)] - 7]ait) + dtF{u{t)) , (22) 



where F{u) is a function of u, and as for the ABBM 
model (HJ 



[F{u) - F{u')]^ =2a\u- 



(23) 



In the limit m -^ the model simplifies to the ABBM 
model considered in the previous section. 

For later use we note that if u(t) is monotonously 
increasing, then the correlator of the random force 
dtF{u{t)) in eq. dH]) is 



dtF{u{t))dt'F{u(t')) = 2au{t)5{t - t') 
= 2<7\u{t)\5{t-t') 



(24) 



Note that ifu{t) > 0, both versions with and without the 
absolute value are equivalent. Under this assumption an 
alternative way to model ([23| is to replace 



(25) 



dtF{u{t)) ^ ^W)m 
= vm\m 



{mm) = 2<7S{t-t'). 



(26) 



We will get back to these formulations shortly. First, let 
us discuss the units. Keeping the same units for time 
and velocity (and space) as given by ([5]) and (O, the in- 
ertial model depends on two dimensionless parameters, 
the driving velocity v := v/Vf^ as before, and the dimen- 
sionless mass 



m/x 



'0 



Indeed two new time scales can be defined: 



(27) 



m/rj 



To = V^ 



m//i (28) 

where t„i is the damping time beyond which damping 
overcomes inertia and tq is the characteristic oscillation 
time of the harmonic oscillator in the absence of damping. 
These time scales are not independent, so there is really 
only one new time scale. The units of acceleration are 

Below we study the model as a function of these two 
parameters m and v. Note that there are various limits 
of interest. We will in particular study the limit of small 
and large m, as well as the limit of large v. The large-m 
limit can be rephrased as the limit where disorder a and 
damping 77 are both small. Note that the weak-disorder 
limit is more general since it is such that only a is small. 

Note a remarkable property of the ABBM model with- 
out inertia: while the space unit ([6]) depends on disor- 
der the characteristic relaxation time ([5]) remains inde- 
pendent of the disorder. Although at small v there is a 
broad distribution of time scales, the characteristic time 
remains t^. In the language of RG it means that the 
friction 77 is not corrected by disorder [2J, [2^]: it is a 
consequence of the Brownian force landscape which can 
be proved using the Middleton theorem. In presence of 
inertia the oscillation and damping time scales tq and 
T„i given in (|28p are the bare ones (i.e. in the absence of 
disorder) and it remains to be understood whether char- 
acteristic oscillation and damping times are affected by 
disorder. 



B. Phenomenology of the inertia model 

The model defined by Eqs. dUll-dMl) is difficuh to ana- 
lyze analytically, since the particle may change its direc- 
tion of motion. Let us therefore start with a numerical 
simulation and some qualitative considerations. We first 
describe how the typical trajectories change as m is in- 
creased, and then we show some numerical results for the 
distribution of instantaneous velocities. 



1. Qualitative features of trajectories 

Let us start with a numerical simulation of Eqs. (HJ)- 
We show in Fig. [T] some examples of typical tra- 
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FIG. 1. Each figure shows two trajectories in the same disorder realization and for the same driving velocity, but different 
values of the mass, as denoted. 



jectories u{t), for different values of the mass and the 
driving velocity in the same realization of disorder, such 
that we can see how the trajectories are correlated with 
the disorder. In the first set we choose a small driving 
velocity v = 0.01 and in the second a larger one v = 0.1, 
so we can see the evolution from the avalanche regime 
(at small v) to the faster driven regime (at larger v). 

For small v we see that upon increasing the mass (up 



to moderate values) : (i) time windows where the particle 
is pinned in a metastable state at position Ui (at zero or 
almost zero velocity) still exist but are shorter since the 
particle oscillates before coming to rest [43|. (ii) Hence 
avalanches (from one metastable state to another one) 
can still be defined, (iii) Due to inertia, the avalanche 
starts more slowly, however the particle overshoots and 
may not settle in the next metastable state (as for m = 0) 



but in one farther away. These metastable states Ui are a 
subset of the metastable states for m = 0. As the mass in- 
creases, more and more of the metastable states at m = 
get eliminated. Thus as the mass increases the smaller 
m = avalanches (i.e. with smaller barriers to the next 
metastable state) are "eaten up" or merge. The larger 
m = avalanches, with larger barriers, remain, although 
the dynamics is quite different. One notes that the fre- 
quency of oscillation increases as the particle settles to 
the new metastable state under the action of damping. 

For larger w as to increases the avalanche structure dis- 
appears and one enters into a regime better described by 
oscillations in the co-moving frame, see Fig. [TJ However 
there remains some correlation with the m = avalanche 
structure, the larger to = avalanches seem to induce the 
largest osciUations (Fig. [T]) . 

While for small mass the motion remains under- 
damped and the time scale remains 0{t^j), in the larger- 
mass regime r^ 3> t^ (to ^ 1 in dimensionless units) 
the motion occurs on larger time scales. Let us describe 
qualitatively the avalanches in that regime. For that it 
is useful to rewrite the equation of motion as: 



-—E = n^vtu 
di ^ 



E 



1 



rju 
V{u) + f^u' 



(29) 
(30) 



where F{u) = —V'{u) and the disorder-potential fluctua- 
tions typically grow as V{u)~V{0) ~ ^v?!"^ . Balancing 
disorder with the quadratic well gives the avalanche cut- 
off size u ^ Sp,. Hence for u <C 5^ the disorder term 
dominates. Consider an avalanche starting at t = 0. 
For small v and ry, as the previous metastable state 
u{t = 0) = Ui becomes unstable, the particle will first 
oscillate with amplitude Au between Ui and the smallest 
root oiV{u) — V{ui) + ^{u'^ — u^), and the total energy E 
can be considered constant during a period J48J . There is 
clearly a distribution of amplitude Aw from the disorder, 
but we can estimate the typical oscillation time Tosc as a 
function of the amplitude Am as 



m'/\Au)^^ya'/^ 



(31) 



by balancing the kinetic energy with the disorder. The 
quadratic well controls the scale of the largest ampli- 
tudes Aw ~ Sfj,, which correspond to a time scale tq. 
The total energy E will decay on much larger time scale 
T,n ^ To , and the particle will settle in one of the avail- 
able metastable states within the range of the first oscil- 
lation. From the above estimate (PT|) one sees that the 
frequency of oscillation will indeed increase as Au de- 
creases to zero. This picture requires that the quadratic 
well has moved by less than S*^ during the avalanche time 
hence v < u,„ := S^JTm (in dimensionless units this is 
V < 1/to). liv > Vm the particle has no time to converge 
to a metastable state and the definition of an avalanche 
becomes less clear. In the regime Vm < v < v^ = S^/tq 
the motion still remains quite correlated to the disorder 



and plateaus would still be visible by averaging over os- 
cillations. Finally for v > vq (in dimensionless units this 
is w > 1/y/rn) multiple oscillations are not visible and 
the trajectory becomes smoother. 

Note that V{u) is a random acceleration process, since 
V"{u) is a white noise. Hence in the large inertial mass 
limit. Aw can be seen as the first return to the origin of 
a random acceleration process. This leads to P{Au) ^ 
^y-5/4 £qj. small Au ^ 5^, the distribution being cutoff 
around 5*^. Similar arguments, although for a slightly 
different model, were made in [2^,|43 using earlier results 

nap. 

It is also interesting to note that in any time window 
[ti,tf] where u(t) > one can parametrize trajectories 
as function of the position u{t) = J dt'u{t') and rewrite 
either model as a stochastic equation for u{u) as 

i!l(lw)H-^(,„, = ,^Q-l)+f». (32, 

In the limit of m — >■ one recovers the standard ABBM 
stochastic equation [171 |l8| and u{u) can be mapped to 
the radial coordinate of a Brownian motion in dimension 
d ^ 1 + v/v^ (see e.g. Section VI B in [52|) leading 
to the distribution of avalanches sizes P{S) ^ S~'^ with 
T = I — Y~ for V < Vfj^, for avalanches much smaller than 
the cutoff, S -^ Sfj^, for which the quadratic well does 
not play a role. In presence of inertia the above equation 
can be used between two zeroes of the velocity and for 
small 77, fjL leads again to an exit time for the random 
acceleration problem. Some considerations about exit 
times are given in App. [H] 



2. Velocity distributions 

Consider on figure [2] the histograms of the distribution 
of the instantaneous velocity ii. We see that for driv- 
ing velocity w = 1, by increasing the mass m from 1/16 
to 16, the distribution becomes more and more symmet- 
ric, and peaked around v — 1. The same happens for 
larger driving velocities, v — 5/2, see figure [3] Note that 
all histograms have a non-vanishing tail for negative ii, 
but that this tail gets smaller when decreasing the mass. 
(Attention: the axis on the figures is shifted from u = 
to the left). By comparing Figs. [2] and [3l we see that 
when increasing the driving velocity v for a fixed mass, 
the probability for negative velocities decreases. This is 
clearly seen on the plots of Fig. [TSl page [231 for mass 
m = 1/4, and v ~ 1/10, v — 1/2, and v ~ 5. The red 
lines (solid on Figs. [21 [3l and dashed on Fig.fTSj) represent 
different approximations to be discussed later. 

This negative tail renders the analytical analysis diffi- 
cult, not per se due to the negative velocities, but since 
the particle moves backwards through the same disor- 
der, thus the disorder it sees becomes correlated in time, 
and the system has memory. Since we possess currently 
no powerful tool to tackle this situation, we will treat 
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FIG. 2. Probability density distribution of the instantaneous 
velocity u for the ABBM model with inertia in the form of a 
histogram for v = 1 and different masses. The red full line 
denotes the probability density obtained in Sec. IIVEI for the 
tree model. Note that the axis is vertical axis is shifted from 
u = to the left. 



FIG. 3. Histogram of the probability density distribution 
for the ABBM model with inertia for v = 5/2 and differ- 
ent masses, in units The red full line denotes the probability 
density obtained in Sec. IIVEI for the tree model. Note that 
the axis is vertical axis is shifted from u = to the left. 



IV. TREE MODEL 



two different local, memory-free variants as explained in 
the introduction: (i) the particle on the tree model that 
can be formulated by using for the random-force correla- 
tor the second line of Eqs. (|24p or ([^5)) : and (ii) the Vu 
model given by the first line of Eqs. p4|) or (|25p . 



There is a strong motivation to consider these two vari- 
ants. Let us look at Figs.[2]and|3l The red full lines rep- 
resent the numerical solution for a particle on the tree. 
We see that when increasing v or decreasing m, the "par- 
ticle on the tree" becomes a good approximation of the 
ABBM model with inertia. The same holds (not shown 
here) for the Vu model. The physical reason is simple: 
At higher and higher driving velocities or smaller and 
smaller mass, the particle will less and less often move 
backward; thus these events will lose their importance 
for P{u). 

We will indeed prove a much stronger statement. Con- 
sider the large-deviation function, defined (supposing the 
limit exists) by 



F(x) :— — lini 



\n[P{xv)] 



(33) 



We will show that F{x) indeed exists for all three models, 
and that for a; > all three large-deviation functions 
coincide. We expect, but cannot prove, that for a: < 
these functions will differ. 



In the next section IIVI we start with the particle on 
the tree, which may be expected to be the most physical 
one. We then continue with the yu model in scction FV Al 
for which we have the most analytical results. 



A. Definition of the model 

In this section we examine the motion of a particle, 
with finite mass m, in a Brownian-correlated disorder 
force. Since the Middleton theorem does not hold, we 
assume that when the particle changes direction, it expe- 
riences a different disorder potential, uncorrelated with 
that experienced previously. The model describes the 
motion on a tree with constraint that the particle always 
chooses a different branch when changing its direction of 
motion. On each branch the disorder satisfies Eq. ©. 
As mentioned in the introduction, the model may be rel- 
evant to describe systems which do not visit the same 
microscopic configuration twice, while the center of mass 
is oscillating back and forth, i.e. systems with large de- 
viations from the Middleton theorem. Note that the tree 
is not defined from the start, but is generated dynami- 
cally, thus there is one tree associated to each trajectory 
or history. This may however be captured by the limit 
of a fixed tree with high branching rate. 

The equation of motion in the laboratory frame is: 



du 

d^ = "' 

m— ^ ^'^[v- u{t)] - rja + dtFt[u], 



(34) 
(35) 



where F is a functional of u{t). Then, the effective dis- 
order correlator becomes 



dtdt>{Ft[u]-Ft.[u]f 



dtFt[u]dt,Ft,[u] 



-Aa\u{t)\5{t-t'). (36) 
2a\u{t)\5{t~t'). (37) 



In the limit of m ^- 0. the model simplifies to the ABBM 
model considered in the previous section. Note that this 
tree model has the pecularity, that the system reaches 



a stationary state without being pinned, even in the ab- 
sence of driving. The reason is that every time it changes 
the direction the disorder is renewed. Only if the system 
starts at a = M = or arrives there "by accident", it 
remains there for ever. 

Now the probability distribution depends on two vari- 
ables, velocity u and acceleration a. The corresponding 
Fokker-Planck equation is a parabolic differential equa- 
tion and has the form 



dP{u,a,t) 
di 



dP 
dii 



da? 

,2 



a\u\P{u, a,i) 



■TT- <, u a H V Flu, a, t) 

oa [ \ m m m 



(38) 



In the following two sections we analyze this equation 
both analytically and numerically. There are several lim- 
its which can be studied analytically, namely small or 
large m at fixed v, and large v at fixed m. 



B. Perturbation expansion in small m 



Equation psp is complicated and an exact analytic so- 
lution is not known. There are many systems [J] where 
the mass term is small and could be treated as pertur- 
bation with respect to the other terms. Therefore, we 
start with perturbation theory in m, i.e. mfj? /rj^ <^ 1 in 
dimensionfull units. 

In dimensionless units the Fokker-Planck equation sim- 
plifies, and effectively depends only on the two parame- 
ters, m and v: 



dP{u, a, t) 



dt 



dP 

ou 


^ da^ \ 


■i\P{u,a,t)\ 
m? J 




-1 
da\ 


( m 

\ m 


a v\ , 
- + -]P{u 
m mj 


,a,i)|. 
(39) 


sted 


in the 


stationary 


situation 



We are in 
dtP{u,a,t) ^0. 

In the limit to — )■ 0^ the acceleration is divergent since 
in the ABBM model the disorder generates a white noise 
with no small-scale cutoff (see Section flV CI below). An- 
alyzing the structure of the Fourier-transformed Fokker- 
Planck equation and the moments that follow from it 
(similarly to Sec. IVPp . we conclude that one has to in- 
troduce a reduced acceleration a = ^/rna in order to be 
able to organize the perturbation theory in to. Then in 
the region |a| > y/m and m > tti > (we will refer to it 
as the region 1) we find: 



P(^)(u,a) = V^e" 



E 



Pn(a,^t) 3±4^ ;, 



nvW2TT 



(40) 



The index (1) denotes that the expression is valid in re- 
gion 1. Here Fn satisfies the recursion for n > (and 
F-i = 0): 

adaFn - iidlFn + u^{v - u)daFn^i + au^duFn-i 
+ -(a^ - {An - l)au)Fn-i = 0. (41) 

We start solving Eq. (|llT) from the smallest n = 0. Then, 
using the solution for F^ we solve the next equation for 
n = 1 and find Fi . The procedure develops further in the 
same way. However, if we want to determine Fi we have 
to solve all the differential equations (|^T|) with n < i + 2. 
The reason is that we a re int erested in a distribution that 
has all finite moments a'^ii^ , and therefore decays faster 
than algebraically for large u and a. This condition has 
to be satisfied for any mass, hence for any order in the 
expansion. Taking this into account when analyzing the 
solution in (i-|-2)nd order allows to discard some solutions 
to the equation appearing in the ith order, i.e. in Fi. The 
details of the calculation and some intermediate results 
arc given in App. [XI 

For brevity we state her only the first three terms; the 
further terms are lengthy, and given in App. [^ 



Fn = l, 



Fi 

F2 



1 

2 V"" 3 

a oa u 

72 ^ 48 

.. 1 



+ C3U , 
1 



(42) 
(43) 



6 



~S -2 

a C3U 



- ~2 • 2/ • 

-a u (u - 



C5U 



48 



'/ [2Au^ - 48uwln(w) - 241;^ 



1~ -3 
-acsu 

36u - 5] . 
(44) 



There still remain undetermined constants C3 and C5. 
They have to be fixed such that J da dii P{u, a) = 1 for 
all TO. This task can not be done now, since region 2 may 
also contribute. It is further complicated by the fact that 
also negative velocities may contribute, and our expan- 
sion does not give a result for those. We discuss this issue 
in the next section, see Eqs. ((9T|) and ((92|) . when com- 
paring the analytical results with the numerical solution 
ofEq. (123). 

Integrating out velocities from P^^'{u,a) one obtains 
the conditional probability distribution of acceleration 
when ii > m. Leading two contributions in the region 
1 are 



P(i'(a)da=da 



2i^ 



2~\ 



\~ar 



K, 



-.(V2I 



\~ar 



3V^r(i;) ■ ■ 
3V2 (a + C3(2i; - 3)) /v |_.„ (^2 |a|) 



2|a|(5-3c3)A'|_„(V2|5|)) 



(45) 



Here Kn{z) is the modified Bessel function of the second 
kind. Similarly, by integrating out a from P^^'{u,a) we 
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obtain the conditional probability distribution of veloci- 
ties in region 1, for \a\ > y/m, of which we state the first 
terms. Higher order terms are discussed in App. El 



P'-^\u) 



C3 



X [(4C5 



T{v) 4r(w) 

1) u - 2u'^ + Auv \n{u) + 2{v 



l)v]. 
(46) 

In the limit of m — > only the first term survives and we 
recover the result given by Eq. (JH) . A note of caution con- 
cerning the distributions (|45p and P5)) is in order: They 
have been obtained by integrating the joint distribution 
(|40|) over all values of a, and over positive ii, while (|40|) 
is valid in a more restricted range. These distributions 
may acquire some correction from other regions, which 
should be small (and maybe even subdominant to the 
correction given above but we cannot prove it). In any 
case, it cannot affect the leading-order result, i.e. the first 
line in P5)) . which is the (normalized) distribution of (re- 
duced) acceleration a in the presence of a small inertia, 
and which is a novel exact result. 

Next we discuss region 2, defined as \a\ < ^Jm and 
< ii < m. If one naively assumes that the results 
obtained in region 1 are valid in region 2, one finds non- 
integrable divergences in P for small velocities u. In order 
to organize the perturbation theory in region 2, we have 
to introduce u = u/m and use a = a/ y/m. Then, 



P(2)(^,a) = /(m,z;)^p(2)( 



a, u)m 



(47) 



Tl = 



where / is an undetermined function. Plugging this equa- 
tion into the Fokker-Planck Eq. (|39p , we obtain the equa- 
tions that determine P„ ' [S^ . Unfortunately, these equa- 
tions are as difficult to solve as the full Fokker-Planck 
equation, and we did not succeeded in solving them 
analytically. However, in App|B] we demonstrate that 
Eq. (|47| has the correct form and that the perturbation 
theory is properly organized. We show that P^'^'{u,a) 
matches at the boundary between regions 1 and 2 the 
distribution function found in region 1. 



C. Large driving velocity v ^ v^ 

A naive argument is that in that limit the noise 
\/\ii\£,{t) can be replaced by ^/v£,{t), hence becomes 
Gaussian. One can then directly solve the Langevin equa- 
tion (in dimensionless units) in frequency space, 

u. = ^^-+'7_f^"^ . (48) 

—raui^ + 1 + zcij 

This leads to a Gaussian distribution for u and a with 

equal-time correlations in the steady state: 

2v 



{u{t)-vf 



~a{ty 



(1 -mcj2)2 -^^2 
2vmuj'^ 



(1 — muj'^Y + ^'^ 



(49) 
(50) 



with u{t)h{t) = and /^ = / §^. We used ^^ = 0. Thus 
we find 



P, 



\,»i(a,?i) w— expf-— [a2 + (u-«)^] j . (51) 



Note that since the dynamics being the same for positive 
ii this large v result should hold for all three models. 

Integrating over the acceleration, this can be identified 
with a Boltzmann distribution 



Pv:»i{u) 



2ttT,, 



(52) 



in the moving frame for the Hamiltonian of a free particle 
of mass m. The effective temperature (in dimensionfull 
units) is 



Tc 



off 



W 



(53) 



Hence at large driving velocity the disorder, the quadratic 
well, and the damp ing act together in the moving frame 
as a thermal noise j54| . 

From the above result we could guess that the proba- 
bility of a negative velocity decays as e""'''-^"''-' at large 
V. We will see below that this is not quite accurate. In- 
deed, we will go beyond the above argument and show 
that there is a large-deviation function which describes 
the deviations from the Gaussian at large driving veloc- 
ity. These deviations appear in the far tails at ju — w| ^ v, 
see Sec. lyTPl 



D. Large mass m S> 1 

At large mass the oscillation time tq = ^/m/ fj, increases 
but becomes much smaller than the damping time r^ = 
m/rj (necessary for damping to overcome inertia) . Hence 
the system is in the underdamped limit and the particle 
oscillates many times before it comes to rest. It can be 
seen by rewriting the Langevin equation in dimensionless 
units: 



du 
dt^ 
da 



,1/4 



m 



a 
m 



(54) 
(55) 



where we have used the reduced acceleration and defined 
the reduced time as t = •Jmil in units of the oscillation 
time. At large m we see that to leading order we have 
a Hamiltonian system with p = a and q = u ~ v and 
H{p,q) = i(p^ + q^), i.e. a harmonic oscillator, weakly 
perturbed by (i) a noise, and (ii) the damping. Although 
these terms are small they will select the steady state as 
we show now. 

For large m the Fokker-Planck equation has a well- 
defined limit if one scales the probability as 

P(u, a, t) — \/rnP{u, a — a^/m, t' — tj \fm) . (56) 
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The scaling function P satisfies 



dt'P = [iii - v)da - adu]P + -^[l + ~ada + D{u)dl]P, 

(57) 
where we have allowed for a general noise function D{u). 
It is equal to D{u) ~ \u\ for the tree model that we study 
now, and D{u) ~ ii for the \/u model studied below. 

This equation is well suited to analyze the time regime 
t ^ To, i.e. t' = 0(1), when the system oscillates. How- 
ever at even larger times it will be damped and will reach 
a steady state. For the latter we are interested in the limit 



lim lim P{u, a,t) . 

771— ^OO i'— fOO 



(58) 



We start by searching for the steady state as a time- 
independent solution in the form 



P{r,e) = ^m-"/2p^(^^, 



(59) 



?i=0 



It is convenient to use the action-angle variables of the 
harmonic oscillator, 



ii ^ V + r sin 9 
a = r cos 9 . 



(60) 
(61) 



In these variables the Fokker-Planck equation becomes 
dgP = ^[l + rcos90 + D(v + rsm9)0'^]P (62) 



O := cos9dr 



sin 9 , 



(63) 



This yields the recursion 

dePo = (64) 

dePn+i = [1 + rcos(6i) O + D(:v + rsin9)0'^]Pn, n>0 

Note that the operators O and D{v + rsin6') commute. 
To leading order the general solution is 



Po(r,0) = Po(r) , 
where Pq it to be normalized as 



/•OO 

/ 27rrPo(^)dr = 1 
Jo 



(65) 
(66) 



/o 
The higher-order terms satisfy 



/ rdr / d6iP„(r,6') ==0, n>\. (67) 

Jo Jo 

The function Pq is selected by the next order equation, 
as we now discuss: 



dePx(T,9,t) 



(68) 



with 



(r, 9) = Po{r) + a{r, 9)P^{r) + b{r, 9)P^'ir) (69) 



1 



air, 9)=r cos^(6l) + -D{v + r sin6') sin^(6') 
r 

b{r, 9) = cos^ {9)D{v + r sin 9) . 



The general solution of (|68p is 



Pi{r,9)^ d9'c^{r,9') + Pi{r,Q) (70) 



r27r 



Now, we observe that one can integrate ((55)) . /^ c\9 and 
obtain a condition which must be satisfied, in order for 
Pi to be meaningful, i.e. a single- valued function in e'^. 



Pi(r,27r)-Pi(r,0) = 0= / A94>{r,9) . (71) 

Jo 

This leads to a condition which determines the steady 
state Po(?') to leading order as 

Po(r) + a{r)Pl,{r) + 6(r)P^'(r) = . (72) 

We have defined 

r 1 f'^'" A9 
a{r)^- + - —D{v + rs\n9)sm^{9) (73) 

2 r Jo 27r 

b{r) ^ — D{v + r sin 9) cos^ (9) . (74) 

Jo 27r 

The first-order correction reads 

P,{r,9) = / d9' [a{r,9')-a{r)]P[,{r) (75) 



+ / d9'[b{r,9')-b{r)]P(,'{r)+Pi{r,0) , 



where Pi(r, 0) is determined from the next-order equa- 
tion using again that P2 should be single valued and the 
normalization condition. Let us now analyze this equa- 
tion for the various models. 



1. D{u) — u model 

Consider now the simpler choice D{u) = ii, which, an- 
ticipating a bit, is the Vu model defined below. In that 
case 



a(r) = - H 



b{r) 



V 



(76) 

(77) 



Solving the differential equation (|72|) we find two solu- 
tions. One is the Gaussian 



^o(r) 



1 

27rw 



-r'^/{2v} 



(78) 



The other one decays as l/r^ at large r and thus cannot 
satisfy the normalization condition (j66p . Going back to 
variables u,a, this identifies with the Gaussian distribu- 
tion ([ST]) also found in the large- velocity limit. 

The analysis can be continued to higher orders. Writ- 



ing 






(79) 
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where the Q„ are polynomials in r, cos 6, and sin 9, which 
satisfy the recursion 

Qo - 1 (80) 

deQn+i = [1 + r cosed + {v + r sme)d^]Qn , n>0. 
We have defined 

6 = e^V(2.)o,-V(2.) ^ eos0 (-'- + dr) —do . 

\ V ) r 

One finds 



Qi 



T cos{9) (r^ cos(26') + r^ - 6v) a^ - 3o 



6u^ 



3w2 



(81) 



The second order Q2 is given in Appendix [Gj The mo- 
ments ii^-a^ of the vu model are computed below in 
Sec. IVDI bv another method and we have checked that 
they agree with the ones obtained here from P to the 
considered order (1/m). 



2. Tree model D{u) = |ii| 
For the tree model one finds 



, . r I— dO 







V 

— h sm ( 
r 



sin2(6') 



h{r) 



1-K 



d0 

2^ 



cos2(6i) 



(82) 
(83) 




FIG. 4. Numerical simulation of the tree model for v = 0.1, 
m = 100, fi — ri = 1, and 5 x lO'' data points. In order to 
eliminate the small asymmetry coming from nonzero velocity, 
we show symmetrized P in A, i.e. [-P(A) + P(— A)]/2. These 
are the data shown (blue dots), together with a la estimate 
of their statistical errors (green, dashed) and the analytical 
curve (|90p (red/solid). Note that there is no adjustable pa- 
rameter. For larger values of |A| (not shown) the statistical 
errors grow and one observe small deviations whose origin 
could be corrections from 1/m or small but finite v. 



We can check that it has zero angular integral and van- 
ishes at all 9 = WK jl. The moments are 



Explicit calculations give 



52 = u2 = i^ + 0(1/771) 

37r^ 



(89) 



37r \r2 



a(r) = -H arcsin 

2 TT r \r 



4|a/1 5", r>v 



«W=2 + 2;' '^"^'^ 


2 + — 






Mr) = -arcsm - + 7— 
TT Vr/ 37r 

V 

b{r) = -, r < V 


y- 


y2 



(84) 



r > V 



(85) 



We see that for large v one recovers the result of the Vu 
model, namely the Gaussian. For v of order one, the solu- 
tion Pq of ([7^ can be computed numerically (see Fig.irU]). 
For V ~ one can obtain an analytical expression using 
a{r) ^ r/2 -f 4/(37r) and h{r) = 2r/(37r): 



97r _3xr 



(86) 



The other solution is excluded, since it decays as l/r^ at 
large r. For the leading correction in 777, at w = we find 

P^{r,9) = P^{r)Q+{r,\e\)sgn{9) , -^ < 9 < t, (87) 



Q^(^'^) = 7^ 4(^-20)(37rr- 
d4 L 



(88) 



-27r 008(6*) (l2r sin[9) + (37rr -j- 4) 008(26*) + Zirr - 20 



whereas they would vanish in the limit of w = for the 
\/ii model, where the Gaussian ()78p is valid at zeroth 
order in I/777. 

The Laplace-transform of Eq. (|86|) gives, to leading 
order in I/777, 



P(A) 



2ir 

d6' 

'0 Jo 

277r3 



drrPo(r)e^" 



(97r2 - 16A2) 



3/2 



(90) 



This function has a branch-cut singularity starting at 
A^ = ±^ = ±2.35619. It is interesting to note that 
from ([7^ and ([M|) one can conclude directly that for 
any u, Po{r) '^ e~^~^ at large r 3> v. Hence the tree 
model has A J = ±^ at tti = cx3 independent of v. Sim- 
ilar branch cut singularities for the ^/u model will be 
discussed in Sec. IV Bl 

In figure 31 we show a numerical simulation of the 
equation of motion given by Eqs. (p4| . (|35|) . and (|36)) 
for V ~ 0.1, 777, = 100, fi ~ rj = 1, and 5 x 10^ data 
points. The agreement of the simulations with the ana- 
lytical result (|M)) is excellent. 
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FIG. 5. (color online) Velocity probability distribution for a 
fixed driving velocity v = 5/2. Different curves correspond to 
different masses, see the explanation in the main text. 



FIG. 6. (color online) Velocity probability distribution for a 
fixed driving velocity v = 1/2. Different curves correspond to 
different masses, see the explanation in the main text. 



Numerical solution of the Fokker-PIanck 
equation 



E. 



In this section we solve Eq. p9l) numerically using the 
discretization scheme proposed by Scharfetter and Gum- 
mel |55| and analyze the probability distribution for dif- 
ferent values of the driving velocity v and mass m. 

The probability distribution of velocities for different 
masses and fixed driving velocity v = 5/2 is shown in 
Fig. [SJ For 771 = there are no negative velocities ii and 
the probability reaches its maximum at a higher value 
than all other curves shown in Fig. [5] The general ten- 
dency is an increase of the probability for negative u when 
increasing the mass as well as a decrease of the maximum 
of the probability. At the same time the maximum gets 
shifted towards higher velocities. We sec that for ?7i = 32 
and m = 128 there is only a small difference in the prob- 
ability distribution, and for large mass the probability 
converges to a master curve at rn = oo, whose behaviour 
we discuss in more detail in the next section. Also note 
that there are two remarkable points where all curves 
intersect, a feature which remains to be understood. 

Fig. |6] shows P{u) for v = 1/2 and different masses. 
We sec that the divergence for m = seems to disappear 
in the presence of a small mass, and that at the same 
time the probability for negative it < becomes finite. 
A small mass changes the positive tail of the distribution 
only slightly. Inertia has the tendency to decrease the 
maximum of P{u) as was noticed in Fig: [5] where v > 1. 
This is also manifest for i) < 1. 

Next wc analyze the probability distribution of accel- 
erations. Fig. [7] shows P{a) for v = 1/2 and different 
masses. The larger the mass is, the more symmetric 
the distribution becomes, in agreement with results from 
Sec. IIVDI Decreasing the mass the maximum of P{a) 
gets lowered, while the distribution broadens and the av- 
erage of a? increases. A similar behavior is observed for 
V >l. 



Apart from similarities between the ABBM model and 
the tree model (discussed in Sec. IIII[) there arc also dif- 
ferences. By looking at Figs. [5] and |31 we see that when 
the probability for negative velocities becomes consider- 
able, a difference between the tree model and the ABBM 
model with inertia becomes visible. The probability dis- 
tribution for the tree model is characterized by a smaller 
peak and larger tails than the ABBM model with inertia. 



F. Comparison of the numerical solution of the 

Fokker-PIanck equation with perturbation expansion 

in m and 1/m 

Here we compare numerical and analytical results from 
the two previous sections focusing on the behaviour at 
small and large m. We start with small m. Wc now 
have to determine the currently undetermined constants 
entering the distribution function, Eqs. (P^ -(P ^ . For 



2.0 
1.5 
I.O 
0.5 
0.0 



FIG. 7. (color online) Probability distribution of accelerations 
for a fixed driving velocity v = 1/2. Different curves corre- 
spond to different masses, see the explanation in the main 
text. 
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sufficiently small m and sufficiently large v, negative ve- 
locities appear with small probability. Then, we can ne- 
glect their contribution to / da du P = 1 as well as the 
contribution from the narrow region 2, since the main 
contribution comes form region 1. We find 



C3 

C5 



=0, 
1 

'4 " 



vil}{v), 



(91) 

(92) 



where ^^{x) is the digamma function. In this region of pa- 
rameters, in App. Klwe state higher-order terms in P{u) 
and additionally we calculate exactly moments charac- 
terizing the distribution function and some other corre- 
lations in Sees. |VDl [VGl and App. |D1 







FIG. 8. (color online) The full green line represents the nu- 
merical solution for the joint probability distribution of ac- 
celeration and velocity, integrated over positive velocities for 
m = 1/8 and v = 5/2. The dashed red line is the result of 
the second order perturbation theory (|45|l . while the dotted 
blue line is the first order of perturbation theory. 

The numerical solution for the joint probability distri- 
bution of acceleration and velocity, integrated over posi- 
tive velocities for m = 1/8 and v = 5/2 is shown in Fig. [8] 
by the full line. The dotted line is the first-order, and the 
dashed line the second-order perturbation theory result 
of Eq. ps)) . We see good agreement between analytical 
and numerical solutions. Also, the agreement increases 
with increasing order of perturbation theory. The small 
remaining difference may come from the approximation 
made when fixing C3 and C5, i.e. due to neglecting the 
contributions from negative velocities and from the re- 
gion 2. 

Let us now compare the distribution of velocities. 
Fig. [3] presents P{u) for m = 1/4 and v = 5/2. The full 
line is the numerical solution, the dotted line is to = 0, 
and the dashed line is the perturbation theory result of 
Eq. (|A6I) . For P{u) in the region 1 there is very good 
agreement between perturbation theory and numerical 
solution. 

Next we compare the perturbation expansion in 1/to, 
studied in Sec. IIVDI with the numerical solution. The 




FIG. 9. (color online) The full green line represents numerical 
solution for velocity distribution for m — 1/4 and v — 5/2. 
Dashed red line is the result of perturbation theory, Eq. (|A6|I . 
Dotted blue line is tti = curve. 



distribution function becomes rotationally invariant in 
the u — a plane around the point (w,0) in the limit 
m — )• 00, see Eq. (p5|) . Therefore, in Fig. [TU] is shown 
P{u+v, a = 0) for fixed mass m = 100 and different driv- 
ing velocities. The definition of P is given by Eq. (j56l) . 
Numerical solutions are shown by full lines, while dotted 
lines denote the zeroth order in the expansion given by 
solutions of the differential Eq. ([7^ for the tree model. 
We see good agreement, especially for larger v. 




FIG. 10. (color online) The full lines represent the numerical 
solution for P{u + v,a = 0) for 771 = 100 and different values 
of V. The dotted lines are the zeroth order of perturbation 
theory, Eq. ^. 



V. Vii MODEL 

A. Strategy 

In this section, we define a model which 

(i) is equivalent to the ABBM model for all trajectories 
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such that the particle only moves forward; 

(ii) is an exact saddle point of the MSR dynamical ac- 
tion hence allows for some analytical results. In 
particular the generating function of the one-time 
velocity distribution is given by 



(93) 



P(A) := c^^M = e^^(^) , 

where Z(X) does not depend on v. 

(iii) on the down-side ii may become complex. 

Ahhough the MSR saddle point method for the ABBM 
model is exact only for m = 0, its extension to m > 
provides a natural approximation of the ABBM model in 
presence of inertia. For instance let us define Papprox('u) 
as the inverse Laplace transform of -P(A), on the real ii 
axis as: 



P(A) 



dU Pa, 



We' 



(94) 



since Z{X)* = Z{X*), the function Papprox(^) is real and 
since Z{0) = (see below) it integrates to unity. As ve- 
locities can become complex, using this function as an ap- 
proximation of a probability distribution (for the ABBM 
model) makes senses if, and only if 

(iv) the function Papprox('") is positive. 

Although not an obvious fact, we have strong numerical 
evidence that this is indeed the case. Finally, we will 
show in Section|VT]that the approximation becomes exact 
for large v. 




FIG. 11. Top: Real (blue) and imaginary (red) part of the 
velocity for one realization of the disorder, and v = 0.5. Bot- 
tom: phase portrait, i.e. trajectories in the complex plane 
u = Ux + iiiy ■ Note the different scale in real and imaginary 
directions. 



B. Definition and basic properties 

We define the model in the laboratory frame, 
du{t) 



At 
da{t) 



a{t), 



(95) 



m—^ = fi'^[v - u{t)] - i^a{t) + ^/{i(f) ^(t) (96) 



{mm)^2<ySit~t') 



(97) 



As long as the velocity is positive, the noise \/u{t) £,{t) 
is real, and the velocity remains real. However when u 
becomes negative, complex velocities arise. This can be 
seen on figure lll[ for one given disorder realization, and 
driving velocity v = 0.5 and m = 1/4. Increasing the 
driving velocity, the trajectories at the bottom of figure 
[TT] get closer to the real axis, and less negative events 
arise. In addition, one sees from ((95|) that near the pos- 
itive real velocity axis a small imaginary part imii = e 
experiences a linear force — (/i^ — /^)^ which on aver- 
age brings it back towards the real axis except for small 
velocities. 



Therefore, this model is expected to be a reasonable 
approximation to the ABBM model with inertia for large 
V and small m, see section lVll Moreover, in the same limit 
of parameters this model works also well for the particle 
on the tree studied in the previous section. It will be dis- 
cussed below in section lV F 6l how good an approximation 
it can provide beyond this range of parameters. 

While the velocities in this model can become com- 
plex, as we will show now, the generating function for 
the velocities is given by a single, v-independent function 
Z(A) 



oA-!i(to) 



,,vZ{\) 



(98) 



The proof is based on the MSR formalism. Replacing 
a(t) in the second line of the above Langevin equation 
by u(i), the expectation ([US]) can be written as [24l - l26{ 



pAit(to) = 



V[u]V[u]e 



r,lpAi'(to)-5[ii,u] 



(99) 
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with the dynamical action 

S[u, u] = 

u{t) [md^u{t) + f]dtu{t) + n\uit) - v)] - au{t)u{tf . 

(100) 

In the square-brackets, we recognize the equation of mo- 
tion without disorder, enforced by the response field u{t). 
The last term resuhs from the average over the disorder, 
using ((57)) . 

Since the action is linear in u, the saddle-point w.r.t. 
this variable is exact, leading to the instanton equation 

md^u{t) - i-idtu{t) + fi^u{t) - auitf = \5{t - to) (101) 

This has to be supplemented with the boundary condi- 
tions u(±oo) = 0. For the case m > studied here it 
turns out to be equivalent to requesting. 



i{t) = y t>to 



(102) 



which implies dtu{t) = for all t > to. [Sg. Supposing 
that (|10ip holds, the only term which survives in ([M|) is 
the term of order v in (jlOOp . which yields ([^S]) with 



z(A) = ^^ys(t). 



(103) 



While analytical results for equations (jlOip and (|103p 
arc difficult to obtain, it is quite easy to solve the instan- 
ton equation poip numerically. We do this now in our 
dimensionlcss units, i.e. setting ?/ = ^^ = a = 1. On 
figure fT2l we show Z{X) for different masses. For real A, 
Z{X) diverges at two singularities, one at A = A+ > 0, 
the other at A = A^T < 0. Their values, as a function 
of m is plotted on figure [T51 Examination of P«(A) for 
complex A show that they are branch-cut singularities. 
The branch-cut singularities determine the behavior of 
the function i'approx(^) defined in ([M|) at large u, up to 
possible power-law pre-exponential factors as 



P. 



approx 



(u) 



e "^<: for u 
QiiXa for ii 



-oo 

-CX3 



(104) 



It is important to note that this asymptotic behavior 
holds at any driving velocity: Were P^„p^o-x^{u) to decay 
faster (e.g. with a larger constant in front of, or a higher 
power of ii in the exponential), then -P(Ac) = c^^^^"'' 
would be finite, which it is not. Similarly, if it would 
decay slower, then -P(Ac) = e"'^^^'^-' would not exist up 
to A = Ac. Note that (|10ip also implies that Z{X) is real 
for real A, at least as long as the solution u{t) decays to 
zero at large times, which is the case for A^T < A < A+ 
(see below). 

As shown in Fig. [131 increasing the mass beyond the 
special value m* ~ 3.95402, the value of A~ decreases 
again, i.e. the tail for negative u becomes again shorter. 
For m — > cx) both A^ and A~ become infinitely large. 




FIG. 12. The generating function Z{X) for masses m — 
(thin, solid, black), m — 1/4 (thick orange), m — 1 (blue, 
dashed), m = 4 (green, dotted), m — 16 (red, dot-dashed, 
thick), m = 64 (black, dotted, thick), m = 256 (green, dotted, 
thick). All functions diverge at A = A^ and A = A^. 



meaning that Z(X) becomes an analytic function for all A. 
This in agrement with the result found in Sec. lIVDl where 
we demonstrate that in the limit m — > oo the distribution 
function becomes Gaussian (f75)) . 

Beyond allowing to define the function Puppro^iu), 
Z{X) contains information on the probability distribution 



P{ux,Uu) of the process u(t) = iix 
plane (see e.g. Fig. [11]) through 



in the complex 



P(A) 



diix 



duj,c^["==+™''lp(u„uj,) . (105) 



It docs not allow however to determine P{ux,Uy) since 
we would need a more general generating function for 
the moments of u* = Ux — iuy A question is what 
information can be extracted, and how does it relate 
to Papprox(u). For iUustration we can consider a sim- 
ple toy example such that the probability factorizes 
P{ux,Uy) = Px{ux)Py{uy) SO that P{X) c&ii. bc writ- 
ten as P{X) ~ Px{X)Py{iX) where Px{X) and Py{X) are 
the Laplace-transforms of Px{u) and Py(u) respectively. 
From our numerical studies we find that P{ux, Uy) is nar- 
rowly distributed in the j/-direction, with a width which 
decreases with increasing v. For large w, P{ux,Uy) will 
converge against Pa.-ppro^{ux)5{uy), at least for iix > 0. 
For finite v consider the simplest example, i.e. a Gaus- 
sian for Py[uy) which leads to 



P,(A)=P,(A)e- 



(106) 



As discussed above, we found that P(A) has two branch- 
cut singularities starting at a finite real A+ and X~ . Since 
the latter control the behavior at large ii and iix, both 
Papprox(w) and Px(ux) have the same behavior at large 
(positive or negative) arguments for this toy example. 
We will not pursue further here the study of P{ux, Uy)- 
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Connection between instanton approach and 
Fokker-Planck approach 



We saw in Section [II] that the ABBM model without 
inertia can be solved cither using the Fokkcr-Planck equa- 
tion, or the MSR method based on the non-linear saddle- 
point equation (the instanton). We recalled their equiv- 
alence based on the method of characteristics. The same 
property holds for the Vw model, with a small caveat. 
The caveat is that the correspondence is simple only 
between the instanton equation and the Fokker-Planck 
equation in its Laplace transformed version Eq. (jlOSp be- 
low. In real space it corresponds formally to the FP equa- 
tion pop where the replacement \u\ — )■ u is made in the 
diffusion term. This would indeed be the FP equation as- 
sociated to the Langevin process ([M]) , (|96p if one forgets 
that ii can become negative, and consequently complex. 
Writing an adequate FP equation for such a process re- 
quires an extension of P{u,a,t) with arguments in the 
complex plane, a route which we do not follow here. 

Instead, consider the evolution equation for the foUow- 



10- 



■ K ■ ■ m^ ■ s ■ ■U's ■ 



-2 



FIG. 13. Top: Location of the branch-cut singularity in Z( A), 
as function of m. Lower curve (red, dotted) is A^, upper 
curve (blue) is Aj. The maximum of A^ is at m, = 3.9540, 
A^ = —3.5124, marked by a dot. Bottom: Blow-up of the 
same curve. 



ing average over trajectories 
P(A, K,i) = e 



Xu{t)-\-K,a{t) 



(107) 



It is in principle an integral / d?i da e'^"~''''°P(u, a, t) over 
the complex u and a plane. The average (|107p satisfies 
the Laplace version of the FP equation. 



dP dP 
dt Ok 



K 

m 



dP 



m 



-vP, 
(108) 



in dimensionless variables. As mentioned above it is for- 
mally the Laplace transform (over the real axis) of the FP 
equation with a ii diffusion term. Eq. (|108p is easily de- 
rived from the Langevin equation (|95p . (|96p by consider- 
ing the variation of the observable in an infinitesimal time 
interval di. The only subtlety arises when expanding the 

^fu{t)da^+o{dt^)to 



variation e > 



1- 



second order, with d^(t)2 = 2adt, using stochastic calcu- 
lus (Ito). This leads to the term {K'^/m^)dxP in pllg)) . 

We now show the connection to the MSR method and 
the instanton equation. The solution of (jlOSp can be 
written in the form 



P(A,K,t) =e"^(^.«^*) , 
where Z is independent of v and satisfies 



dZ dZ 



dt Ok 



dZ 






d\ 



(109) 



= — (110) 
m 



To solve this equation we again apply the method of char- 
acteristics, and reduce the partial differential equation to 
a family of ordinary differential equations. 



k{t) = -\{t) + 



K{t) 



m = ^ - <'^" 



m^ 



dZ(i) = ^dt. 
m 



(111) 
(112) 
(113) 



We have defined Z{t) = Z{X{t), ^(i), t), as well as k{t) 
dK{t)/dt and X{t) = dX{t)/dt. 
Let us now introduce 



u(t) = K{t)/m 



(114) 



Eliminating X{t) in (jllip and pi2p we find that it satis- 
fies 



mu{t) - u{t) + u{t) - u\t) = 



(115) 



If we impose u{t) ~ and u{t) = for t > to and as- 
sume A(to) = A and K(to) — 0, this equation becomes 
the instanton equation (|10ip in dimensionless units. To 
obtain the stationary solution Z{X,k) of piOp we can 
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solve (|115p on the interval t G] — oo,t*], with boundary 
conditions at i = i* : 



u{t*) 
ti{t*) 



A K 

— + ^7 



(116) 
(117) 



. 2 
<M > 



0.80 
0.75 
0.70 
0.65 



= 1/2 



-4-2 2 4 



Logjm 



<M > 

9.1 f 

9.0 

8.9 



8.7 
8.6 



=5/2 



-4-2 2 4 



Logjm 



and u(— oo) = 0. u(~oo) ~ 0. Then we compute 



Z{t*) 



u{t)dt ;= Z{X,k) 



(118) 



FIG. 14. Average value of u^ for the ABBM model with iner- 
tia as a function of mass m is shown for two driving velocities 
V = 1/2 and v = 5/2. Dashed line represents the value v{v+l) 
while the points are the numerical results. 



which is precisely Z{X, k) if expressed as a function of 
the boundary condition. Since it does not depend on t* , 
we have found the stationary solution. 

Hence wc have shown, via Eq. ()108p . that the observ- 
able (|107p can be obtained from the solution of the in- 
stanton equation, although with a slightly more general 
boundary condition than in (|10ip . This is because we 
now want the joint distribution of velocity and accelera- 
tion (at a given time t*). We can rewrite the observable 
pU7)) as 



function for the 



model: 



P(A,K,t) = el <it[^s{t-t')-K.S'it~t' 



)]u{t) 



(119) 



Performing the same mani pula tions using the dynamical 
MSR action as in Refs. |24| - (26| as sketched above in 
Section IV B[ we arrive at the instanton equation with a 
source on the right-hand-side XS{t—t*) — K6'{t — t*) which 
is equivalent to the boundary conditions (|116p . (|117p . 

To summarize, for the V^ model the equation (|108|) 
describes the time evolution of the observable (|107|) un- 
der the Langevin equation (|95p - (|97|) . even though u and 
a take values in the complex plane. Solving it is equiv- 
alent to solving the instanton equation as a function of 
A and k that determine its boundary conditions (|116p 
and (|117p . Neither equation admits solutions in closed 
analytical form for generic values of m, but each has his 
advantages for numerical or perturbative studies. For 
instance, from Eq. (jlOSp wc can easily obtain the mo- 
ments for this model, as we now show. Also, note that 
the Laplace transform of the Fokker-Planck equation for 
the tree model can also be studied although it is more 
involved, see App. O 



D. Moments of the distribution function 

In the stationary case, taking the derivatives (?"9™ 
of Eq. pOSp where m and n are positive integers and 
afterwards setting k and A to zero, one obtains a set of 
linear equations for 9J^9™P(A, k)|(o,o)- By solving them 
and using 9;'a^'P(A, k)|(o,o) = u^^a" and P(0,0) = 1 



i g g - 



v,o„i 



u"a^ = 0, 

V 



2v 



m{m -f 2) ' 
Su^ 6(5m + 3)v 



m? m2(m + 2)(4m + 3) 



0, 



2v 



m m{m + 2) ' 
2v^ 



6v 



m{m + 2) r?T,(m -|- 2)(4m -f 3) ' 



ii^a*^ = v(v + 1), 



0, 



1 



{Am V + Am V 



2„,2 



m(m-|-2)(4m + 3) 
-I- llmv^ + 27mw^ -I- lOmw + 6i;^ + 18w^ + 12w), 
^5— ?r 4ti 



ii'^n^ 



u-'a 



m + 2 
0. 



+ w^ + 3w^ 



(120) 
(121) 

(122) 

(123) 

(124) 
(125) 

(126) 

(127) 

(128) 
(129) 



(130) 
(131) 

(132) 



one finds [57'| the moments characterizing the distribution 



For brevity we stated only the first few moments, but 
the procedure can be easily extended to higher moments. 
Some of them are given in App. [D] 

While these results are exact for the vu model we can 
compare them with the ABBM model with inertia and 
for the particle on the tree model. In Fig. [14] we show 
(m^) as a function of mass for two different driving ve- 
locities for the ABBM model with inertia. For v = 1/2 
we see disagreement between the above obtained v{v + l) 
value and the numerical result for m = 1, while for the 
larger velocity the deviation appears at larger values of 
the mass. Similar behavior is observed for the tree model. 

In the App. [D] the moments of the velocity are com- 
puted using an iterative solution of the instanton equa- 
tion instead of Fokker-Planck equation. 
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E. Perturbation at small m for P{u) 

111 this section we obtain the perturbative expansion in 
small m of the generating function Z{X) = Z{\^ n ~ Q) 
for the steady state velocity distribution function of the 
yii model. To this aim we will solve perturbatively in m 
the non linear instanton equation (|115|) and use (|118|) to 
obtain Z{\) from the solution. 

The solution u{t) of (|115p at ttt, = is given in (|16p . It 
is non-zero for i < 0, zero for i > with a jump at < = 0. 
In presence of inertia tti > we are now looking for a 
solution non-zero for < < 0, zero for t > 0, but which 
vanishes at f = 0, i.e. m(0) = since we have set k = 0, 
and with u(0) — —\/m. For small m it is then clear that 
there is a boundary layer for t — 0{m) < where the 
function varies rapidly from values of order to values of 
order A. In addition there is a bulk region t = 0(1) < 0. 
We now study both regions and how they match. 

As will become clear below, in the boundary layer re- 
gion < — i < 771 we can write 



l{t) = Y, "^"/n - 



(133) 



with m(0) = and u{0) = —X/m. Inserting into the 
instanton equation (|115|) we then find the following re- 
cursion: 

/o(<) = A(l-e*), (134) 

k 

fk+i{t) - h+i{t) + fk{t) - Y, Mt)fk-dt) = 0, (135) 

1=0 

for fc > 0. Note that the term fo{t) already exists in 
absence of disorder and describes the rounding of the 
disorder- free response function by the mass. We assume 
/„(0) = and /„(0) = forn > 0. Then we find 



/iW=^|4-2t(A-l) 



5A 



A^ 



o2t 



e*A 2A-2-t(2A- 1) 



(136) 



Higher-order terms can be found, but for brevity are not 
state here. In the App. |E] we give /2. In general /„ 
is a sum of exponentials and a polynomial, of the form 
Ut) = YZ\ e'=*C("'(t, A) + ELo At\X){-tf further 
discussed in App. |e1 

In the region — t = 0(1) ^ m the solution can be 
written as 



i(t) = ^m"y„(t) 



(137) 



n=0 



where ?/„ satisfy the following differential equations 



- yo + yo - 2/5 = 0> 
yk-i{t) -yk{t) + yk{t) 



e=o 



(138) 
yi{t)yk-i{t) = 0. (139) 



The last line is for k > 0. The boundary condition is 
y„(— oo) = 0. Eqs. (|138p and (|139p can be solved recur- 
sively. This gives for each yn{t) a first-order differential 
equation, for which we need to fix yn(0) by the condi- 
tion that solutions in different regions (|133ll37p match 
at — t small but —t/m large. We discuss this in detail in 
App. El Here we state the solution 



yoit) = 
yiit) = 



A 



(140) 
+ 5A - 4(-l + A) In (e^* + A - e~*A)] . (141) 



A + (1-A)e-*' 
e-*A 



Additionally, in App. |E] we give y2{t)- Now we can de- 
termine 



Z(A,0)=^m"Z„(A) 



(142) 



given by Eq. pTS)) with t* = 0. We find Zq = 

f^^dtyoit) and Z„ = /!^ dt ELi e'*ci"-^)(i. A) + 

/_ dty„(t). Hence there is a contribution from both 

regions. The contribution of the Aj^ terms is already 
taken into account through the contribution of the y„ (t) 
at small t. Calculating the integrals we obtain: 



Zo = - In (1 - A) 
Zi 



2 Vl-A 



21n(l- A) 
A 



Z. = ^(^Mi±^)^M±M+i,(i_,). 



24(A-1) 



(143) 
(144) 

(145) 



NowP(A,0) ^ e^ 



ju2(A,o) g_j^(^ from that follows the 



perturbative expansion of the moments it*-' in agreement 
with the results in Sec. IV Dl Assuming that for small 
mass we can ignore contribution in P coming from com- 
plex velocities we perform the inverse Fourier transform 
of P and find the distribution function 



= (") 



-diu) 



me '^u" ■^ 



T{v) ' ' 2T{v) 

X (— 2uv In (w) -I- 2uviIj{v) + ii^ — v^ + v) 6{u), 
+ 0{m^) (146) 

where 0{x) is the Heaviside theta function, and V'(^) is 
the digamma function. The result to this order is thus a 
probability normalized to unity (to this order) for v > 1. 
The correction 0(m?') is given in appendix[^by Eq. (|A6p 
and is a bona-fide distribution for v > 2. It also gives 
nonzero contribution for u > only and it is normalized 
to zero. Note that there is no a priori reason why at a 
fixed TO > the inverse Laplace transform of Z{X) should 
exist, and furthermore with support on m > since the 
present model leads to complex velocities. However, it 
appears that (i) for large enough driving velocity v > n 
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the perturbativG result can be trusted up to order 0{m"'), 
i.e. the negative velocity events are sufficiently rare not 
to spoil the result to that order; (ii) the perturbative 
result seems to be correct for any v provided u > m. 
These findings are consistent with the analysis for the 
tree model, for which we found the same result (|46p in 
the region u ^ m although with unknown constants C3 
and C5 while for the vw model they are fixed by Eqs. 
(lOTMl) . 



F. Critical case rn = 1/4 

Here we analyze the specific case of ?ti = 1/4, and 
work in our dimensionless units. For the system with- 
out disorder this is the critical case which separates the 
over-damped to under-damped regime where the particle 
starts to oscillate. It is also the mass which brings the 
particle back to the origin of the parabola in the shortest 
time. In presence of disorder there is a priori nothing spe- 
cial to this case. Since, however, perturbative formulas 
simplify for m = 1/4 our aim in this Section is to obtain 
to high numerical accuracy the behaviour of Z(X) and of 
Papprox{u)- While the results are specific to m = 1/4, 
they are representative of other masses. 

Since neither Z(X) nor Pappiox(w) can be obtained ex- 
actly, we consider several schemes to compute them to 
excellent precision (i) a perturbative expansion in A; (ii) 
an analysis of Z close to the branch-cuts, and (iii) the be- 
havior of Z for A — >■ 00. Each of these schemes gives an 
estimates of Papprox(w), (the latter two are are denoted 
Pbc{u) and Pasymp(w) below) depending on which region 
of the complex plane dominates in the Laplace inversion. 
They arc valid for different values of v and u, as discussed 
below. In addition we use numerical solution of the in- 
stanton equation. At the end we compare our results to 
the numerical simulations on the ABBM model. 




FIG. 15. Z{\) in the critical case m = 1/4, together with its 
Taylor expansion (|149|) . 
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+ 



7657689375 71693475840000 
289088854220889511A3 
" 5357932381952640000 
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+0 (Ai3) (149) 

The branch cuts described in Section fV Bl are at 



A+ = 1.10647 



Ar = -8.58563 



(150) 



1. Perturbative expansion for Z{X) 



2. Perturbative expansion for the critical instanton and X^ 



The response function is 



R{t) = 4te-2*e(t) , 



^ + 1 ) R{t) = 5(t) (147) 



The instanton equation (|10ip becomes 



The value of the branch cuts can either be obtained 
from the numerical solution of the instanton equation, 
or via the following observation: The critical instanton. 



ut{t) 



(t) converges for t 



—00 rapidly against 
1, see figure [TBI Thus while Z(A^) = 00, the solution of 
the instanton equation is still defined at this point. This 
allows for a series expansion, making the ansatz 



1 - ^ ) ut-ui = XS{t) 



(148) 



The boundary conditions are given by Eq. (J102I) . Eq. 
p48p can be solved iteratively in A, as shown in App. |D] 
Integrating over time yields a perturbative expansion of 



-{t + t*) = l-c° 



)i=2 



bne" 



(151) 



The parameters a > and 6,1 have to be determined, 
such that (|15ip satisfies the instanton equation (|101|) . 
and finally t* is chosen such that u'^(t*) = 0. To find 
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FIG. 16. The function ut is increasing with increasing A, as 
seen by numericaUy integrating the instanton solution: A — 
1/2 (green, dashed, smallest curve), A = 1 (blue, dashed), 
A = 1.106 (black, dashed), A = A+ (red/blue, thick). The 
last curve shows superimposed the numerical solution of the 
instanton equation, and of the expansion (|152|) . 



the parameters, expand ()15ip inserted into the instanton 
equation (|101|) in a series in y := c"'. At first order 
in y, the condition is 4 + 4a — a^ = 0, leading to a = 
2(1 + v2). Then solve order by order in y, determining 
the 6„. Finally, find t* , such that u+{t = 0) = 0. (This 
is necessary since the higher orders shift the origin.) The 
result for the critical instanton is 

u+{t) = 1 - 1.08835y + 0.0935856y2 - 0.00550998?;3 
+0.000284501y^ - 0.00001369y'^ 



+6.30567 X IQ-'^y^ 



2.81909 X IQ-^y^ 



-1.23335 X 10~V - 5.30802 x 10""^^ 
-2.25516 X 10-i2yio _ 9.48217 x 10""y" 
-3.95296 X IQ-^^y^^ 



0{y'-') 



(152) 



This determines A^ = 1.10647. Of course, all coefficients 
could be given explicitly. Also note that already for n = 2 
one gets A^ within an error of 0.01. Each further order 
in n improves the accuracy of Ajt by about one order 
of magnitude. Thus one could use this algorithm as an 
efficient estimator for A^, as a function of m. 
A similar solution could be constructed for Ar . 



3. A good approximation for Z{X) from the branch cuts 

There is an astonishingly good approximation for Z{X), 
which is dominated by the branch cuts: 

Zbc(A) = -1.20711 In ( (^-K)i\->^t) \ (153) 



It is obtained by starting from the masslcss case, 
Zjn^Q^X) = — ln(l — A/A|j") with A+|„i=o — 1; a-nd adding 
a similar term — ln(l — A/A~) for the negative tail. The 



prcfactor was then determined by asking that the highest 
known moment, in (|149p of order A^^ be given exactly. 
Then 

Z{X) - Zbc(A) == 0.0496465A- 0.00117287A2 

-0.0000986924A3 - 2.769730882007382 x lO'^A^ 

-f-2.48119 X lO^^A^ - 5.57763 x 10~"A^ 

-2.96177 X lO^^A" - 1.82040 x lO^^A*^ 

-1.04881 X lO^^A^ - 5.56368 x 10"^°A^° 

-2.25950 X 10"i°A" + O (A^^) (154) 

What is remarkable about this approximation, besides its 
accuracy in reproducing higher moments, is that both ex- 
pected terms, ln(l — A/A+) and ln(l — A/A^t) appear with 
the same amplitude. Performing the inverse-Laplace 
transform for f > one finds, defining v — 1.20711w: 



= -i«(AJ+A+) 



ftc(u) 



^-'^>^^r{jf^y 



xK~, 



V^r(w) 

i[A+ - A-] ) e{u > 0) 



5-iu(AJ+A+)(_^-^+yr. 



A^-A+ 



xK._ 



^u[K 



-Xt])eiu<0). (155) 



K is the Bessel- iiT function. The function Pbc {u) decays 
as e""^": " for u to oo, and as e~''*<: " for ii to — oo. When 
compared to the numerical inverse-Laplace transform, we 
find that this is a good approximation for all but small 
ii. 

In the limit of ?j = 0+, we know (see [2^ for more 
details) that one can expand P(A) = e"^'-^'^ = l+vZ(X) + 
0{v^), and that upon Laplace inversion one obtains 

P{u) = (1 - vp')S{u) + vp'P{u) -f 0(w2) . (156) 

The first term represents events when the particle is 
pinned, and the second one yields the velocity distribu- 
tion in an avalanche P{u) via Z{X) = wp'Jdu(e'''" — 
l)P{u) where vp' is the probability that t ~ to belongs 
to an avalanche. From the form (|153p for Zbc(A) we get 



p'Pbc(u) = 1.20711 



-,-ua: 



^— uA„ 



-e{u>o) + 



(u < 0) 



(157) 

Note that for m > it is very similar to the result for 
the ABBM model with m = [111 up to the different 
value for Aji" . More interestingly, it also gives a non-trivial 
prediction for the tail on the negative- velocity side in the 
avalanche regime. 



4- Large \-behavior 



In order to obtain the small-w behavior of Pj 
one has to analyze Z(X) for |A| — >■ cx). For real A this 



approx V^ji 
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is impossible since branch-cuts appear at A = A^. As 
we show in appendix |Fl the behavior for A — >■ ±ioo on 
the imaginary axis can be calculated analytically, and is 
given by 



^asymp 



(A) 



4.76006 (|A|to)^/^ . 



(158) 



The contribution of expression (|158p to Pasymp(w) at fi- 
nite driving velocity v is 



Pasympiu) = / COs(A|u|)e 

Jo T^ 

■ cos(A|u|)e 
tV3^exp( °-"°387W. 



dA 

7^ 



vZ{\i) 



|u|V3 

1.43193m2/^uAi 



4.6923|up/3Ai 



0.715967m4/3w2 
., /^0.715967m4/3i;2 



(159) 



Ai is the Airy function, and the result (|159p a positive 
function, peaked around zero, even for relatively large 
driving velocities, see figure [TT] Its value at w = is 



asynip 



(0) = 



0.181215 



(160) 



Note that while p59p is relevant at finite w, it does not 
contribute to the large-deviation function discussed be- 
low. We will discuss below its domain of applicability. 




FIG. 17. The contribution of the large-A behavior (|158p 
Pasymp(M) at w = 2, m = 1/4, as given in Eq. (|159p . 



to 



First the perturbative expansion in A works well for 
large v. This will be quantified in section IVTl where we 
discuss the large-deviation function. 

Second, the branch-cut approximation (jl55p is a rea- 
sonable approximation to -Papprox(^) for all v and all ii, 
and becomes a good approximation in the tails. The lat- 
ter can be expected since it uses the knowledge about the 
branch-cut singularities at A^ and A~. 

Third, the approximation Pasymp(u), given in p59p . 
It is a reasonable approximation at small u, as long as 
i; is small enough. For v = 1/2 e.g. it predicts, with 
relative precision of 10^'^, the value of [Papprox(0''') -I- 



-'approxt.U )\/ ^ 

from 1. 68172 foru = 



where we note that Pappiox(u) jumps at 0, 



to 2.42249 at u = 0+. The worse 
disagreement at ii = 0.1 is probably due to problems in 
the numerical inverse Laplace-transform (P(A) oscillates 
strongly on the chosen contour). For w = 5, it does not 
work. 



5. From Z{\) to Pa, 



.(«) 



The instanton equation can be solved numerically for 
any complex A away from the branch cuts on the real axis 
at A > A+ or A < A^. This yields Z{\) for complex A. We 
have performed the numerical inverse Laplace transform 
for Z{\). A convenient choice of the contour is A = 
a(l ± i), a > for u > and A = — a(l ± i), a > 
for u < 0. This gives an excellent numerical accuracy, 
except when ii and v are both very small (< 0.1). 

We can now compare with the asymptotic estimates 
of Papprox('") discussed above (see Figure [TH]). Why the 
different approximations work or fail, can be traced to 
an analysis of the inverse Laplace transform. Depending 
on ii and w, it is dominated by one of the three special 
points: A = for the perturbative approximation, valid 
for V — >■ oo, A = A^ for the tails, and A — >■ ±ioo, for small 
M, as long as v is not large enough s.t. A = dominates. 



6. Tlie ^/u-model as an approximation for ABBM 

On Figure [T8l we show the data obtained for the prob- 
ability distribution P{u) from a numerical simulation of 
the ABBM model at m = 1/4 (green shaded area), com- 
pared to results obtained in this section. The driving ve- 
locities are v = 1/10, v = 1/2 and v = 5 (from top to bot- 
tom). First compare with the result for Papprox('"), ob- 
tained by numerically inverse-Laplace transforming P{X) 
(thick red dashed line). The agreement is excellent for 
V — 5, and reasonable for v ~ 1/2 (it does not give well 
the peak for ii close to 0, but quite well the tails, even 
on the negative side. For w = 0.1 only the tails are rea- 
sonably well approximated. (We cannot however exclude 
numerical problems in the inverse-Laplace transform for 
it < 0.01.) 
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m 




m = 6/25 for which there exist analytic solutions in 
closed form. These are known in the context of the 
Fisher-Kolmogorov equation and Abel equations, which 
is intimately linked to our problem as we now explain. 



Fisher-Kolmogorov equation 



Eq. (|115p is a second order differential equation and it 
has many solutions that are characterized by different set 
of values (kqj Aq) define by Eqs. (|116lll7p where t* is an 
arbitrary time. One solution for m = 6/25 reads as|58j: 



u{t) 



1 



(l + i?e-5*/6)2' 



(161) 



where D is an arbitrary parameter. Here we used that 
traveling wave solutions of the Fisher-Kolmogorov equa- 
tion satisfy Eq. (|115|) , and that for a special value of the 
speed of the propagating wave it has the analytic solution 
Eq. (|16ip that corresponds to the special value m = 6/25 



m our case. 



FIG. 18. Comparison of Papprox(M) obtained by numerically 
inverse-Laplace transforming e"'^*-'^' (thick red dashed lines) 
with a simulation of the ABBM model (green shaded his- 
togram plot), the approximation (jlSSp (black dotted lines), 
and Pasymp('J) giveu by . The mass is m, = 1/4, and the driv- '^I'^j 



For D > it follows that < ii < 1 and using Eq. (jllSp 
we get 

Z{X{k),k) = -- (\n (1 - ^25k/6) + y^25K/6] , (162) 

we determine below which values A takes as a function 
of n. For D < solution u{t) (given by Eq. ()161|) ) takes 
some value that is greater than 1 at two different times. If 

one chooses t* = -6\n((l + y/6/{25K)) /{-D)) /5 then 

Z{X{k), k) = -| fin (1 + ^25k/6) - ^25^/6') . (163) 

Here u takes all values greater than zero. If one takes 

t* = -6\n((l-y/6/25^] {-D))/5, then Z{X{k),k) = 

oo because the integration in Eq. (|118p is over the di- 
vergence that happens at t = 61n(|£'|)/5. Taking into 
account these results one finds the following correlation 
functions 



ing velocities are from top to bottom: 
V = 5. 



1/10, V = 1/2, and 



gX(K)u+Ka ^ pVZ{\(K),K)^ (164) 

'5k/2 4- 50k3/2/63/2 Z? > and < k < 6/25 
5k/2 - 50^3/2/63/2 !)< and K > 0. 

(165) 



G. Exact results for m = 6/25 



While Eq. (|115p cannot be solved in closed form for 
generic m, and one had to resort to expansions or nu- 
merical solutions, there exists a magic value of the mass 



In Eq. (|164p . the first and the second line of Eq. (|165p 
correspond to Z{k) given by Eq. (|162p and (|163p . respec- 
tively. Hence the Laplace transform of the joint distri- 
bution of velocities and acceleration is known exactly on 
the curve A — A(k), which, upon expanding in k can be 
translated into non-trivial relations between moments of 
this distribution. 
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2. Abel differential equation of the second kind 



Next, we find more general solution of Eq. (|115p than 
in the previous section. Then using it we obtain an exact 
result for Z(A, 0) for some range of A values. Addition- 
ally, we obtain an analytic expression for A^ that tell us 
about behavior of the tails of P{u) for positive velocities, 
sec Eq. (fTOT]) . 

Introducing A = X/rn — K/rn? and B = k/tti, the set 
of Eqs. (|111I112I113P can be rewritten as: 

dA B^-B-A , ^ 

AB mA 

dZ _ B 

Eq. (|166p is the so-called Abel differential equation of the 
second kind. Its parametric solution reads as |59| 

1 



Bc{s) = -^E\s,C)s, 



(168) 



Ac{s) = -— S2(s,C) (Vl + s3£;(.s,C) + 2s) , (169) 



E{s,C)^ / dr(l + r3)-i/2^C, 
Jo 



(170) 



where C is an arbitrary constant. We see that i?c(0) = 
while ^^(O) = — C'^5/24 can be arbitrary, meaning 
that we have found the family of solutions determined 
by parameter C . Then 



Z{\c{s),Kc{s)) 



Ac{s) ds 



(171) 

where the constant is to be determined by knowing that 
Z(0,0) = 0. 

For example, if we want to calculate 

exp (-C3u/20) = cxp (wZ(-CV20, 0)), (172) 

we need to find first 



^'-^^'•»>-£/' 



Bc{s) dBcis) 



Aci-s) ds 



= -ol d-^ 



=Eis,C), (173) 



where tco is defined as E{tco,C) = 0. Wc find analytic 
expression for 

^(-^C-^, 0) =^ {t'co + 1) 3^2 (^l, 1, 1; I 2; t'co + 1 



50F V6/ V3, 

1 1 3 
3'2' 2' 



x\Aco + l2Fi(^,^;^;4o + l 
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C^co ^Fi 



12 5 3 
2' 3' 3'"*^" 



-^H'-"(S) 



(174) 



Z(A,0) 



15 



10 



(167) FIG. 19. Numerical result for Z{X, 0) is shown for m = 6/25. 



Here 2^1 is the hypergeometric function and 3F2 is gen- 
eralized hypergeometric function, tco is implicitly given 
by E{tco, C) = that could be rewritten as 



Eitco,0) = -C. 



(175) 



We find that E{s,C) = C + S2F1 (l/3, l/2;4/3; -s^). 
The function E{s,0) becomes complex for s < — 1. For 
s > — 1 it is increasing function of s, and satisfies: 



r(i)r(^) , , r(i)r(|) 



(176) 



We see that for Amin < A = — C^/20 < Amax, where 



1 
20 



2V^ 



Amax - 2Q 



mm 



-0.137843, (177) 



1.10274, (178) 



there exist a solution of Ec^. (|175p and one can find tco- 
Since at s = tco both A(s) and k{s) arc zero, this con- 
straint on A might mean that a solution of Eq. (J115p 
with boundary conditions K(t = 0) = and k{t = 0) = A 
reaches {k{t),K{t)) = (0,0), i.e. it is convergent meaning 
that at large times it goes to zero. By considering nu- 
merically this instanton solution of Eq. (|115p with given 
boundary conditions, one indeed sees that it starts to 
diverge for this value of A and A+ = Amax- However, 
for Amin the instanton solution is still convergent. The 
numerics gives A~ = —8.8219 < A,„in = —0.1378. In 
Fig. [19] we show the numerical result for Z(A, 0). We 
conclude that the parametric solution (|168ll69ll70p does 
not "catch" whole K{t) dependence on time with given 
boundary conditions K{t = 0) = and k{t = 0) = A for 

A < Aniin- 

Next, wc solve E{tco: C) = pcrturbatively in C and 
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find 



tco = -c ^ 



C^ 



Ci 



Z(A,0) = A+ — + 



cr_ 

112 ^ 1792 

75C19 

20873216 793182208 '^ 

A2 25A^^ 125A4 

- + 
2 

375A5 



3C 



13 



93184 



84 616 

115625A6 



2548 



1039584 



+ 0{X'). 



(179) 



(180) 



Now, expanding Eq. (|172p in A = — C"^/20 one can find 
the moments u'"" in agreement with Sec. IV Dj and App.ID] 



VI. LARGE-DEVIATION FUNCTION 

A. Definitions and numerical determination 

It is suggestive that we "do a good job for large driv- 
ing velocities" , since when the particle does not move 
backward, all three models considered here are indistin- 
guishable. We actually show below that the so-called 
large-deviation function coincides for all three models for 
X > 0. The large-deviation function is defined as the 
leading behaviour for large driving velocity v of the dis- 
tribution of instantaneous velocity P{u) as follows: 



F,ix) 



In [P{xv)] 



F{x) := lim F^{x). 
Analogously, one defines, if that limit exists, 



^.(A) 
Z(A) 



In e^"^ 



(181) 
(182) 



(183) 
(184) 



lim Zy{X) . 

i;— >-oo 

If the limits exist, then for large v the Laplace-transform 



gA« ^ g.vZ^{X) 



dxe-^'-^-t^^e^'"" 



(185) 



can be approximated by its saddle point, and Z{X) and 
F{x) are related via a Legendre transform: 



Z{X) + F{x) =Xx 
AZ(A)^.. 



(186) 
(187) 

(188) 



It is assumed that F" [x) > 0. Note that Z{X) is easier 
to measure numerically than F(x), since the former does 
not need binning. 

Let us now review our numerical results and how they 
are consistent with the following scenario: 

(i) Zy(X) becomes u-independant at large v for each 
model in some range of A. 



(ii) the asymptotic curves coincide for A > A*, 

Z(A) = Z^{X) = Zabbm(A) = Zt,cc(A) , (189) 

where 

Z'{X*)=0 (190) 

i.e. the minimum of Z{X) is at A*. This corresponds to 
the point of zero velocity since Eq. (|188|) implies a;(A*) = 
0. Another way to state Eq. (|189l) is to say that for .t > 



Fix) 



(191) 



-^\Ai(^) = -PABBM(a;) = -Ftrco(a;) • 

We assumed that Z'{X) > for A > A*. A simple argu- 
ment, which shows that (|189p and (|19ip holds, is given 
below. For A < A* the Z{X) for each model is dominated 
by negative velocities. Since the models differ for these 
velocities, there is no reason why their Z{X) should be 
the same, and consequently -^(2;) for a; < is expected 
to be different for the different models. 

Our numerical data for Zy (A) are presented on the left 
of Fig. [5D] for m = 1/4. First of all, we have checked 
through large-scale simulations (5 x 10^ data points) that 
for the aAx- model Zy{X) = Z{X) is velocity- independent, 
and given by the numerical solution of the instanton 
equation (jlOip . using (|103|) (thick orange line). Only 
data for w = 0.1 are presented on the plot (dark grey- 
blue line, with error bars given by the dashed lines of the 
same color). Relative errors are < 10% for —4.2 < A < 1. 

We have then checked that 2'„(A) for the tree-model 
converges, for w — >■ 00, towards the numerical solution 
of the instanton equation (jl03p from above, while the 
ABBM model converges from below. (This fact is con- 
sistent with Figs. [2] and [31 where one observes that the 
probability distribution for the tree model has larger tails 
than for the ABBM model with inertia.) On Fig.[20l data 
for V = 0.5 and v = 2 are shown for the tree model, and 
for w = 5 for the ABBM model. Simulations for other 
velocities (not shown) confirm this picture. Note, how- 
ever, that we have found convergence of the simulations 
of Z„(A) for both ABBM and the tree model to Z{X) of 
the %Ai-modcl only in some domain A > A* and certainly 
not for A < A*, corresponding to negative velocity, con- 
sistent with the above scenario. For m — 1/4 we find 
A* = -4.26953. 

Once the Z{X) curves have been measured for each 
model, one can obtain the corresponding large-deviation 
functions via a numerical Legendre transform. On fig- 
ure [ini we have plotted the large-deviation function, 
Legendrc-transform of Z{X) for the yu model. Its min- 
imum is at .T = 1, and for large positive x it grows like 
F{x) ~ A+x -I- const, with A+ given by equation (|150p . 
For large negative x, the growth is F(x) ~ X'^ x + const 
with A^ also given by equation (|150p . These slopes are 
indicated by the dashed curves in figure [20] (right). 



B. Convergence of the large-deviation functions 

We now show the main result of this section, namely 
that the three models have the same large-deviation tunc- 
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FIG. 20. Left: Different curves for Z{\), m = 1/4: The numerical solution of the instanton equation (|10ip is used to obtain Z 
(see p03p ) denoted as vu inst and shown by a thick orange line. Numerical simulation of this \Ai- model for driving velocity 
11 = 0.1 (grey-blue line), with an estimate of the numerical error bars, (one-cr error, dashed curves of the same color). Within 
error-bars, the simulation has converged to the orange curve Vii inst; errors are small for —4 < A < 1 . The simulations for 
the tree model at u = 2 (dark green) and v = 0.5 (bright green), both with error bars (dashed) show clear deviations from 
instanton solution at small v, but get closer to the latter for larger v. The last curve is a simulation of the ABBM model at 
V = 5 (red). It has also converged to the orange curve within error bars. This is not the case for smaller driving velocities (not 
shown). All simulations are for 5 x 10^ data points, apart from the one for ABBM, which has 93.000 data points. Due to the 
small number, the error bars in that case are underestimated. 

Right: The large-deviation function F{x), obtained by performing a numerical Legendre transform of Z{X) given by the 
numerical solution of the instanton equation (fTOT|) . and (fTCTSl) . We note F(0) = F(3.47268) = 1.30459. The thick red part is the 
domain for which F{x) < F(0), which can be obtained with arbitrary precision from a simulation with negligible probability of 
negative velocity as explained in the text, hence must coincide for the three models. The dashed curves show the asymptotic 
behavior XcX + const. 



tion in the region < a; < x+i 

F{x) = F^{x) = Ftrcc(a;) = -FABBM(a;) 

for all < X- < a;+ , (192) 

F{x+) = F(0) (193) 

The argument is simple and in fact much more general: 
Two models, which have the exact same dynamics for 
positive velocities, have the same large-deviation function 
in that interval. The idea is the following. Consider a 
simulation with a set of Ny = e'^'" data points. If d > 
F(0) there will be negative velocities in a typical such 
set (with probability one at large v). But if rf < F{Q) 
there will be none, again with probability one at large 
V. Consider any x such that F{x) < F{0). We can then 
measure the value of F{x) with arbitrary accuracy (as v 
becomes large) if we use e''" data points with F{x) < d < 
F(0) (in fact d = i^(x)+ is sufhcient). As stated above, 
this set almost surely does not contain negative velocities. 
Since the dynamics for the three models exactly coincides 



for trajectories with positive velocities, this shows the 
above property (|192p . provided at least for one of the 
models the large-deviation function exists. However the 
latter is true for the yu-model, since there Z{X) is v- 
independent, which completes the argument. 

The argument is based on considerations of under- 
sampling of rare events, and is reminiscent of similar con- 
siderations used for the multi-fractal spectrum of wave- 
functions, e.g. when comparing the size-dependence- 
exponents of participation ratios for a typical sample or 
disorder-averaged ones I60,l61[. Here, the additional in- 
put is the identification of the dynamics for positive ve- 
locities. Note the restriction that x < x-i-, which means 
that rare events with positive velocities but as rare as 
the negative velocity can not be controlled neither. In 
the present case we do not see a reason why the various 
functions F{x) would not coincide for x > x+. The re- 
striction comes from the generality of the argument: One 
could for instance imagine a dynamics such that the par- 
ticle jumps discontinuously from large positive velocities 
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to negative ones. The above estimates are made more 
precise in appendix HI 



C. Large deviation function in perturbation theory 
for small m 



Expression obtained from previous results 



We note that F{x) = F{x,y = 0), which presumably 
holds to all orders. 

One can check that convexity F"{x) > holds for x > 
77i/(l + m) for the two first orders in m. This indicates 
that the expansion is valid at small m only for x 3> 0{m). 
We will see below, on the example of to = 1/4, that 
indeed the small-m expansion is accurate for large enough 
x and breaks down for small x, see figure [H] 



In this section we compute the large-deviation function 
for the Vu model in a pcrturbative expansion at small 
m. For a; > 0, as argued above, it gives the result for all 
three models. A further discussion of this equivalence is 
given in App. [Cl 

We slightly generalize the discussion of the previous 
Section by considering the large-deviation function for 
both velocity and acceleration defined as 



F{x,y) 



lim 

^— >-oo 



In [P{xv,yv)] 



(194) 



where x = u/v and y — a/v. The connection to the 
generating function introduced in (|1071 1109|) is through 
a Legendre transform in both variables, namely 



Z(A, n) = max{— i^(x, y) + \x + kj/} 



^,v 



(195) 



To compute F(x,y) from Z(A,k), one looks for (xo,2/o) 
that satisfy dxF\xo,y„ = A and dyF\xo,y„ = k, to get 



F{xo,yo) = -Z{X, k) + Xxo + kj/o, 



(196) 



These formula assume somehow the convexity of F(x, y) 
and Z{\, k) which we did not attempt to prove but seems 
to hold. 

Now let us use our results from Sec lIVBl to obtain 
F[x, y) perturbatively in m. To this aim we use Eq. (jl94p 
with P{u,a) given by Eqs. (gO]) and (l42])-(l44]) where C3 
and C5 are determined by Eqs. (|9T|) and (|92|) . and u is 
restricted to positive values. We find 



v' 



F{x,y) = [^+x-\n{x)-l 
' 2x 



my^ 



6a;2 



f y' 



V48a;3 4x2 ^ 4a; 2 2x ^^''' 



+ m 



3/2 



f 



f 



^\ 



240x4 36x3 132,2 ; 



+ 0(m2), 
(197) 



where y = \fmy. 

Similarly we obtain the large-deviation function for the 
velocity only, defined in Eq. p82p . Legendre transform- 
ing the results of Sec. IV El for Z(\) for small to, using 
(IT86)) , we get 

( X 1 
F(x) = X — ln(.T) — l + TTi ln(x) 

1 5x 3 , , X 1 

T2^-l2 + 4^+^"("')~4 



^(to^). 



2. Equation from Fokker-Planck 

From the Fokker-Planck equation it is possible to ob- 
tain a differential equation for the large-deviation func- 
tion. Inserting the form P{u, a) = e""^"^^^^'^^"-' in the 
Fokker-Planck equation ([M)) we find that F{x,y) satis- 
fies, to dominant order at large u, 

rn^yd^F - to(x + y- l)dyF + D{x){dyFf = . (199) 

Here D{x) ~ x for the ^/ii model and D{x) ~ \x\ for the 
tree model. If we study this equation for x > 0, the two 
equations are the same. While the equation for the tree 
model can be studied for all x, the meaning of the one 
for the vw model for x < requires further analysis due 
to possible complex velocities. 

We now use equation (|199p . and the emerging struc- 
ture of the above pcrturbative results (|197p and (|198p , to 
construct the expansion in m to higher orders. One way 
to analyze (jl99p is to perform a Taylor expansion around 
2/ = 0, 



F(x,y)=Fo(x) + ^F„(x)y" 

n>2 



(200) 



where we assume F{x) = F{x,y = 0). Then the 
Fn{x) obey the recursion relations Fi{x) = 0, ^^2(2^) = 
F^(x)/[2(x - 1)], and for fc > 2 



Fk+i{x) 



to(x- l)(fc + l) 



kFk{x) + ^/^iFl_^{x) 



fc-i 



+x^{j) + l){k - p + l)Fp+i{x)Fk-p+i{x) . 

(201) 



If we assume that the structure of (|197p holds to higher 
orders, in particular that there are no poles in l/(x— 1), 
we find that the conditions for their cancelations order 
by order in m give enough conditions to determine Fq (x) 
entirely. The form which we find by inspection is 



Fk{x) = Y.' 



n" + 



p=n+k-l h 

k-2 X — > O-r, 



E"'np 
tP 



=/c-l 



XP 



(202) 



(19(3) j^Qp k > Q. One checks that the previous result (|197p 
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satisfies equation (|199p . One finally obtains 

F{x) = Fo{x) = {x- \n{x) - 1) 

1 / 1 , , , 
+ —m Ix 2 In(a-) 



+ m ln(a:;) — 



5x'^ -9x + l 1 
12;^?2 4 



a; 107x^ - 240x + 22 +5 



21n(x) 



53 



I 144x3 '"'"' 72 

'669 - xix{U809x^ - 38610x + 2090) + 2905) 



8640x4 

^, , , 3895\ ^, 5, 
51n(a;)-^^j+0(m'^). 



(203) 



We have fixed the integration constant by requiring that 
F{1) = 0, a consequence oi ii = v. Note that the above 
result satisfies F'{1) ~ and F"{1) = 1, which is consis- 
tent with our analysis in Section IIVC| namely that the 
bulk of the distribution is the Gaussian ([ST|) in the scaling 
region |u — u| ^ ^/v. The first corrections to the Gaussian 
arise from F"'{1) = — 4/(2 + m) and may already be visi- 
ble in the tail |u— w| ^ v"^^^. This can be compared to the 
perturbative expansion of F{x) around x ~ 1, given in 
Eq. (|207p below. The complete large-deviation function 
describes the far tails |m — w| ~ v. Another interesting 
feature is that one obtains an expansion for X'^{m) as 



A+(m) 



Fix) 



lim 

a;— i-oo X 

"^^y ~ ~[2^ 144 



14809m'' 
8640 



(204) 
+ 0{m^). 



Finally note that for large x, y the large deviation func- 
tion takes the form: 



F{x,y) Kxg{y/x) 



(205) 



and that an ordinary differential equation can be written 
from dUSl) for g[z). 



D. Large deviation function for large m and 
matching of small and large m 

In App. C it is explained how to calculate the pertur- 
bative solution of Z{\) in powers of A. Examining the 
results in powers of A keeping the complete m depen- 
dence of the coefficients shows that it actually turns into 
a \aige-m expansion. Indeed one finds 



z(A)^A+^+^^+ ^;(^-+^) 

^ ' 2 3to + 6 16m2 -|- 44m 

2A5(103?7i2 + 198m-h72) 
5(m + 2)(m -f- 6)(4m + 3)(9m -I- 4) 
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Af^(6957 



4396m3 + 7666m2 + 4284m + 720) 



3(m 4- 2)2(m + 6)(4m + 3)(9m + 4)(1677i + 5) 
O(A^), (206) 



and each two orders more in A come with an additional 
factor of \/m at large m. We have obtained two more 
orders, which are not displayed here due to their length. 
Legendre transforming yields 

2 6m + 6 

(12 +16771- 57^2) (x- 1)4 
^ 4(m + 2)2(4m-|-3) 

4(61m4 + 420m3 - 338m2 - 540m - 144) (x - 1)^ 
^ 5(m + 2)3(m + 6)(4m + 3)(977i + 4) 

{x-lf 

^ 6(m + 2)4(m + 6)(4m + 3)2(9m + 4)(16m -f 5) ^ 

X [(17280 -t- 143136m + 386448m2 -I- 239312m3 

- 488936m4 - 556346m^ - 3195m^ -t- 5240m^) 
+ 0(x-l)^ (207) 

which is actually an expansion in x— 1 at fixed tti, in other 
words deviations from the Gaussian solution of Sec. IIVCI 
at large velocity. Again, we have obtained two more or- 
ders, which are not displayed here due to their length. 
Since it yields the derivatives i^^")(l) for any ttt,, one can 
check that at small m, they match the result obtained 
from (pg5| above. 

From the above expansion one can now obtain a good 
approximation for -F'(O), which gives an estimate of the 
probability 



P' 



-vFifi) 



(208) 



for negative velocities at large w, and improves on the 
estimate of Sec. IIVGI 

As a test, we can compare the small- 777, expansion (j203p 
and the large-777, expansion (|207p with the numerical so- 
lution of the instanton equation, followed by a Legendre 
transform (in the form of a parametric representation). 
The result for 777 ~ 1/4 is shown on figure W\\ One sees 
that the small- 777 expansion works well for large x, but 
breaks down for .x — > 0, while the large- 771 expansion con- 
verges for X = 0, but may have a finite radius of conver- 
gence in X— 1. Taking both expansions together, we have 
an analytical approximation for the range < x < x_|_ 
drawn in red in Fig. [20] in its right part of Sec. IVI Al 
where the large-deviation functions for the three models 
have been argued to coincide. 



E. Large-deviation function for m = 6/25 

Here we determine exactly, in a parametric form, the 
large-deviation function for the special value of the mass 
777 = 6/25 and for a certain range of values for x, using 
the results obtained in Sec. IVG21 

One finds F(x > 0), using Eq. (fTM)) . where Z{\ = 
-CV2O) is determined by Eq. (fT7i)) . In Eq. (|W)l . x 



29 




FIG. 21. F{x) for 771 — 1/4 as obtained by the numerical 
solution of the instanton equation (dashed line). Large- tti 
expansion (|207|) to order 1/m^, thus two more orders in (l — x) 
as given in Eq. (|207|) (blue solid line); this approximation 
substantially deviates for x > 2. Small m expansion (|203|) . 
which works well at large x, but breaks down at about x = 1/4 
(red solid line). 
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(209) 



where ico is given by Eq. (|175|) . This expression is real 
as long as tco > — 1. 



VII. THERMAL AND QUANTUM 
FLUCTUATIONS 



As long as we neglect backward trajectories, it is pos- 
sible to include thermal and quantum fluctuations. Here 
we obtain some new results within that approximation. 
We then discuss its expected range of validity. 



A. Classical systems: Thermal fluctuations 

For a classical system in presence of a thermal noise, we 
can generalize the equation of motion in the laboratory 
frame ^, ^ as 

mu{t) + 7ju{t)=fi^[v~u{t)]+dtF{uit))+^T{t) ■ (210) 

To remain slightly more general, and cover the case of a 
colored noise, we define the correlations of the noise as 



The standard thermal white noise is 



B{t) = Brit) = 2r]Td{t) . 



(212) 



Hence when we can neglect the backward motion, and for 
F(u) a Brownian landscape, the problem reduces again 
to a Langevin equation for the velocity, 

mu{t) + i]u{t) = /i^[w - u{t)] + Viii,{t) 4- irit) (213) 

in presence of both a thermal noise and a noise generated 
by the quenched disorder. It still describes the equation 
of motion for the center-of-mass velocity of an elastic 
manifold moving in a Brownian force landscape in pres- 
ence of a thermal noise, under the same (but more strin- 
gent) approximation of neglecting any backward move- 
ment of the interface (i.e. even of a piece of it). Indeed 
for forward motion, the general argument given in Refs. 
m Hi] and [H] (section IV.A) still applies. 

Also note that the Langevin equation (|213p defines a 
noisy version of the y/u model, as a legitimate model 
provided one accepts complex velocities which arise from 
backward trajectories. 

It turns out that ()213|) can be solved exactly for an 
arbitrary noise correlator B{t). This is remarkable, since 
except when ^{t) is a white noise itself (see below) no 
simple Fokker-Planck version seems to exist. However, 
within the MSR formalism the problem is much simpler. 
The Laplace transform of the velocity distribution can 
be written as in Eq. (|99|) . where the dynamical action 
S\u,u\ contains the additional term 



u{t)u{t')dtdt'B{t-t') (214) 

dtu{t)dt'u{t')B{t - t') . (215) 



The derivation follows the same steps as in Section IV Bl 
The action is still linear in m, hence the instanton equa- 
tion pOip is unchanged and inserting this equation in the 
action gives 



^g/ A!i(to)) ^ gVZW + h fit' dMt)d,,u{t')B(t-t') (^21^) 



using that u{t) vanishes at t 



Z{\) 
gives 



' j^u{t) as in p03)l . 



= ±00, with the same 
For a thermal noise it 



^g/Au(to)) ^gfZ(A)+,,T/jatfi(t)]^ 



(217) 



{iTmT{t'))^B{t-t') 



(211) 



Note that a double average is performed over thermal 
fluctuations and disorder realizations [63 ■ Note also that 
these results can be obtained from the general expression 
given in Ref. [2q| (see Eq. (6) there), valid for an arbitrary 
forward driving and a forward motion w{t), by substitut- 
ing w{t) — > vt + £t(0//^^ ^i^d performing the Gaussian 
average over ^t |63| . 

The calculation of the contribution of the noise requires 
a small-time cutoff. For standard thermal white noise 
with no intrinsic cutoff, (j217p can be evaluated only for 
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a non-zero inertial mass, which provides the small-time 
cutoff. Its evaluation turns out to be remarkably sim- 
ple, as we now show: Multiplying the instanton equation 
(jlOip by dfU, it can be rewritten, for t < to, as 



-dt 



W{u) 



2 



W{u{t)) 



'r^[dtu{t)? (218) 
(219) 



This can be interpreted as a classical particle of position 
variable u{t), undergoing (backward in time) a damped 
motion in a potential W{u). Integrating over time we 
find 



dt{dtu) 



2 

2m 



{dtu{t)? + W{u{t)) 



lO 



(220) 



using that u{to) = and u'{to) = A/m, as well as that 
u(t) decays to zero at i = — oo. In other words, the 
total dissipated energy is the total initial energy since 
the fictitious particle settles to rest &i t = — oo. Hence, 
we find the exact result for any m and d. 



^g/Au(«o)) ^e"2(A)-l-, 



(221) 



The new term corresponds to the thermal broadening 
(i.e. u — > u -|- 5u) of the T = velocity distribution by a 
Gaussian of variance {Sii^Y = ^- ^^ '-^'^ ^^ interpreted 
as the system reaching kinetic-energy equipartition, as 
in the system without disorder. Keeping only the terms 
of 0{\) and 0{\^) in (|22ip leads to the large-w Gaus- 
sian distribution (|46p with Toff -^ TcS + T, and indicates 
that at large v negative velocities arc negligible when 
T + Toff ^ mv'^ . More precisely, a new large-deviation 
function at large v can be defined if one scales T ^ v, 
i.e in the high-temperature, high-driving-velocity limit. 
The Legendre transform of Z{X) +9X^/2 yields to Fg{x). 
Defining = T /(rnvvV) we have 



^-vFi){x=^) 



P{u) 

F'{x)=F'g[x + eF'{x)). 



(222) 
(223) 



For the case m ~ Q (ABBM model without inertia) we 
find using Fg{\) — 



Fe{x)^{x-l) 



-0 



+ 9+ sjm -H (x - 



In 



^[x-9+^J\9+{x 



(224) 



The second solution that follows from p23|) is not ap- 
plicable since it does not give Gaussian distribution ex- 
pected at large w, Sec. I IV CI This formula should, fol- 
lowing the general argument given in Section I VII also 
hold in the positive velocity region a; > for the orig- 
inal ABBM model (without inertia) and with tempera- 
ture, even though it cannot be solved exactly because of 



possible backward motion due to thermal fluctuations. 
Extension of (|224[) in presence of inertia can be studied 
using the array of methods introduced in this paper, but 
we refrain from doing so here. 

The result (|22ip is remarkable since the effect of ther- 
mal fluctuations is exactly Gaussian, despite the presence 
of quenched randomness which is highly non-linear. This 
property can be generalized and traced to the fact that 
the noise-dissipation satisfies the fluctuation dissipation 
relation (FDR), i.e. for an equilibrium thermal bath. To 
see that consider a slightly more general bath and re- 
sponse function, in frequency space. 



R ^{uj) = —muP' + /i^ + ri{uj)iuj , 



(225) 



where rj{uj) is an even function in u. The classical FDR 
reads 



2Tt]{oj) = B{u) . 



(226) 



Since the response function is changed, the instanton 
equation becomes different. It involves the transpose of 
the inverse response and reads 



{mdf + ii)u - 



au 



*ij{Lo)iLou{uj) = . (227) 



Again, multiplying with dtu{t) and integrating over time 
yields 

= / 7/(w)a;^u(a;)i7(-a;) . (228) 

2to J,., 



Hence 

1 /■ A^ 

dtu(t)dt'u(t')B(t -t')^—- (229) 

2m 

if the FDR holds. Eq. (PTCl) again implies Eq. ((^17)) . 

Formula (|216p also allows to compute e-*'"+'"', and the 
joint distribution of velocity and acceleration, replacing 
Z(A) by Z{X, k) and u{t) by the solution of the instan- 
ton equation with boundary conditions (|116lll7p . The 
formula (|217p however cannot be used, as it contains a 
divergent integral since u(t) has a jump of k/tti at to- 
This is because the distribution of acceleration is not well 
defined; unless we add an intrinsic small-time cutoff to 
the thermal bath, i.e. the mass is not sufficient to act 
as a cutoff. It can be seen within the large-w analysis of 
Section I IV CI since (a^) is defined by the same integral 
with an additional lu'^ in the numerator. 

A non-Gaussian contribution to the velocity distri- 
bution can arise from a Gaussian noise if it is a non- 
equilibrium one, such as, for instance colored noise with a 
frequency-independent constant-friction dissipation term 
rj. We mention here one simple example, when f (t) is a 
white noise, {i{t)i{t)) = 2DS{t - t'), i.e. ^(i) is a Brow- 
nian. Then it is easy to see that for m = one recovers 
the result for the ABBM model without inertia ^ for 
P{u) with the replacement 



P^,^d{u) = AP^+D^D=o{u + D) 



(230) 
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This is the only case which is amenable to a Fokker- 
Planck equation (with u — > u + D in the diffusion term) . 
Not of course, that negative velocities are simply ne- 
glected, hence the normalization factor A. 



B. Quantum system 

The extension to the quantum system in presence of 
a bath is straightforward. It can be done by generaliz- 
ing the MSR methods of Ref. [13, ^ to the Keldysh 
path integral. Let us write the Keldysh action for a 
quantum particle in a Brownian-force landscape. The 
Keldysh path-integral over the fields u{t),u{t) between 
an initial and a final time involves e~^ with Keldysh 
action Sk = S'^ + S]^ and 






W Jt 
1 



iu{t)B{t-t')iu{t') 
w 



(231) 
(232) 



^k = '^\Y. ^("W + 4^W - ^W)' (233) 

^i .- LI 



£ = ±1 



iu{t)[?{u{t) -w{t)) 



(234) 



The path integral can be expressed alternatively in terms 
of the upper and lower Keldysh fields u^{t) = u{t) ± 
|u(i). The form of the functions GQ^{t) and B{t) de- 
pends on the details of the bath. One convenient choice 
is 



Sk + SI^ I lutimd'i + r^dt - tJ.'')ut 



(235) 



-- / iutB{t-t')iut'+ i 
^ Jw Jt 



lUt^JL Wt 



with 



B{uj) = rihiocoth{j3hLo/2) , 



(236) 



which for 77 > represents an Ohmic bath. A realistic 
bath has a large- frequency cutoff ujc- 

After averaging over disorder the system is described 
by the same Keldysh action with S]^ replaced by 



ee'R\ut - Uf + -{eut - e'lif)] . 



"'« e,e' = ±l 

(237) 

This is because the Keldysh path integral is normalized 
to unity. 

In the classical limit fi, — >■ 0, one recovers the classical 
MSR functional with the thermal white noise B{uj) ~ 
2i]T and the usual disorder part 



SI 



classical 



1 



iutiuf^iuxt - Uxt') 



(238) 



A(u) = —R"(u) is the correlator of the pinning force. 
Note that iii is sometimes denoted by u in the MSR for- 
malism. 

Now we can treat the case of a Brownian force land- 
scape choosing 

i?(u) = ii-(0)-A(0)y+aM^. (239) 

This corresponds to A(u) = A(0)— cru. Inserting this into 
Eq. (|237|) . we obtain a complicated expression. However, 
if we make the replacement 



sgn{u\t) ~ u'' {t')) ^ &gn{t - t') 



(240) 



for all four couples (e, e') = (±1,±1), then it simplifies 
into 
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1 



UfUr 



A(0) - f7(ut - utOsign(t - t') 



The observable we are computing is the following average 
over the Keldysh action 



P[\]^{eM^'^^)s^ 



(241) 



The study of more general observable is left for the future. 
To recover the velocity theory we define 



lUt 



-dfUt , 



(242) 



and consider ut vanishing at t = ±cx). It is very similar 
to what is done in [2Q|, and to which we refer for details. 
It yields, after integration by parts 



Sk = -a 



UtUfUt 



(243) 



Sk + Si-= / ut{md^ + Vdt - fi^)ut 



(244) 



dtUtB{t - t')dt'Ut' - / iutWt 



Integration over iit leads to the instanton equation ()10ip , 
and inserting its solution into the action, we find again 
for w{t) = vt the same result (|216|) . Note however that 
for any non zero h the bath cutoff time is needed to get 
a finite result - the mass only cutoff leads to a logarith- 
mic divergence when inserting into (|216p . Since they re- 
quire the corresponding dissipation related by FRD, we 
leave explicit calculations to future studies. To summa- 
rize however, one can say that everything works as if the 
quantum system is described by a semi-classical equation 
of motion with the noise correlator (|236|) . 

Of course, there are two crucial ingredients here: 
(i) the Brownian disorder-force landscape; (ii) the ap- 
proximation (|240p . For e' = e it amounts to neglecting 
any trajectory with backward motion. We see however 
that for e' = — e the approximation cannot be correct at 
short time differences t — t' , even for forward-only tra- 
jectories. This should be valid if this time scale is much 
smaller than the other ones considered here. A more de- 
tailed discussion of the validity of this approach is left for 
the future. 
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VIII. CONCLUSIONS AND DISCUSSION 

In this paper we studied an extension of the ABBM 
model including inertia, i.e. the motion of a particle of 
inertial mass m driven at externally imposed average ve- 
locity V > in a ID Brownian random-force landscape in 
presence of damping. Its main interest, besides modeling 
a particle, is that in some cases it describes the center-of- 
mass dynamics of an interface in a Brownian correlated 
disorder. When all the segments of the interface move 
forward it is certainly true, and in more general situa- 
tions it remains to be understood. For any m this model 
can also be derived for the center of mass of a manifold 
with inertia in a short-range disorder potential, in the 
limit of a fully connected model (infinite-ranged interac- 
tions). In that sense it is a mean-field approximation. 
Whether that property also extends to finite range in- 
teractions in high enough internal space dimension d, as 
is the case for m = for d > dc, remains to be under- 
stood. Our aim here was to calculate exactly, or using 
approximate methods, the distribution of the instanta- 
neous velocity, and of the acceleration for this particle 
model. 

We started by recalling the 771 = limit which is ex- 
actly solvable (standard ABBM model). It is charac- 
terized by a relaxation time scale t^, a spatial scale S*^ 
and a velocity scale w^ = S^/r^, see Sec. |lll For v < w^ 
the motion proceeds via avalanches, while for v > v^ 
the motion becomes smoother. In presence of inertia 
TO > 0, the ABBM model cannot easily be solved be- 
cause backward motion occurs and induces memory ef- 
fects. Numerical simulations and qualitative arguments 
showed new characteristic time scales for oscillations tq, 
and damping t™, and new velocity scales vq = S^/tq and 
Vm = S^/ttti fSec lIIl] ). An avalanche regime survives for 
1; < niin(i;m, u^) similar to the one for to = except that 
the smallest avalanches have merged into bigger ones. As 
the mass increases, overshoots and oscillations become 
more pronounced before the damping allows relaxation 
into a metastable state. As v increases further the mo- 
tion becomes smoother, but oscillations persist. As a 
general rule the inertia tends to make the motion less 
jerky and to smoothen the abrupt jumps of a centcr-of- 
mass position in time. At the same time, the distribution 
of velocities evolves from being strictly positive but with 
a divergent limit P(u = 0"'') = cx3 for ttt, = and u < 1, 
corresponding to intermittent motion where the particle 
is part of the time at rest, to developing a finite weight for 
negative velocities and a finite P{u = 0), corresponding 
to oscillatory motion. 

To make quantitative progress we introduced two vari- 
ants of the model which share the exact same dynamics 
with the ABBM model for all forward trajectories, and 
are analytically more tractable. The analysis of the first, 
the tree model, is based on a Fokker-Planck equation. 
The study of the second one, the y/u model, is based 
on a saddle-point equation of the dynamical action. For 
both of them we calculated the joint distribution of ac- 



celeration and velocity pcrturbativcly in small and large 
mass. The y/ii model could also be solved exactly for 
a magic value of the mass in terms of hypergeometric 
functions, and was studied with very high precision for 
other values of the mass. From these variant models we 
obtained two sets of results for the ABBM model: 

(i) at large driving velocity: Since by increasing v the 
probability distribution for negative velocities decreases, 
all the considered models become more and more similar. 
The bulk of the velocity distribution, i.e. for velocities 
|u — iij = 0(w), then tends to the same Gaussian. To 
characterize more accurately the tails we defined for each 
model the large-deviation function which describes the 
rare events when the instantaneous velocity ii deviates 
from the average v. We proved that these large-deviation 
functions become identical for positive u at large v for 
the three models, and obtained analytical expressions at 
small and large to, and for the magic value of the mass. 

(ii) at any driving velocity: We compared the three 
models and discuss differences and similarities. Although 
agreement is not exact anymore, some features of the 
ABBM model are, in some cases, quite well reproduced. 

Finally we showed how thermal and quantum fluctua- 
tions can also be treated within the approximation of ne- 
glecting backward trajectories. For thermal fluctuations 
it is expected to be a reasonable approximation at flxed 
V only for small T or for any T at large v. For quantum 
fluctuations, the discussion is more subtle but basically 
it should hold within a semi-classical approximation. 

In conclusion this paper proposes a first step to the de- 
scription of classical and quantum avalanches of pinned 
elastic systems in presence of inertia, and the effect of 
driving. A more elaborate theory should incorporate bar- 
rier crossing by thermal and quantum fluctuations, and 
a treatment of memory and oscillation effects. However 
we believe that we have introduced a useful framework. 
For instance a key observable is Z{\) ~ ifiie^" which 
in the V" model is independent of f , and is well charac- 
terized by branch-cut singularities describing the tails of 
the velocity distribution. At the same time, as w — )- 0+ it 
describes the avalanche dynamics. Hence, numerical or 
experimental determination of this quantity in realistic 
systems could help developing further understanding. 

It would be interesting to extend the current analysis 
to different time-dependent and space-dependent driv- 
ing, as well as to analyze the spatial correlations of the 
probability distribution taking into account the spatial 
extension of the interface and not just its centcr-of-mass 
position. Additionally, due to retardation effects appear- 
ing in soft magnets [J] there is a need to generalize the 
current approaches and to use a more general equation of 
motion, with a more general response function. Finally 
in quantum systems where velocity translates into cur- 
rent, developing a more general connection with the full 
counting-statistics problem would be very interesting. 
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Appendix A: Perturbation theory at small m for 
P{u) for the tree model 

In this appendix we give details on the perturbation 
theory discussed in the Sec. IIVBI Wc start with the dif- 
ferential equation (|^T|) with n = 0. It has the solution 



Fo{a,u) 



l-^'''^''^''''[vm 



+ C2(u), (Al) 



where crfi denotes the imaginary error function crfi(z) = 
erf(iz)/i. In order to have properly defined moments 
a^ii^ , the distribution function has to decay exponentially 
fast at large a and ii. This implies that ci(w) = 0. Next, 



Fi(a, u) = Ci{u) 



1 

+ 2 



-—a^C2{u) + auc2{u) — 2au^C2{u) 



(A2) 



F2{a, it) ~ -a^ii^ [— wc2(w) + uc2{u) — uc2{u)] 

+ 3ra'^C2(u) 4- -j^a^u [8uc2{u) - 5c2{u)] 



a^Ci{u) 



+ -rci^u^ [2u^C2{u) + VC2{u) — UC2{u)\ 



^Vilc,iu)cr& (^ 



V2v^ 



-ail [2iic^{ii) — 5c4(-u)] 4- ce{ii), 



(A3) 



where 2F2{a;b] z) is a generalized hypergeometric func- 
tion. The first line of Eq. (|A3[) has to vanish due to its 
large-a behavior, which gives us a differential equation 
for C2('u). Solving it, we obtain 



C2iu) = C3(-1)'T(1 - w, -u) + a 



(A4) 



Analyzing it, we conclude that C3 = 0. Then, finally we 
obtain that FQ{a,ii) = C4. Since the distribution has to 



be normalized by J dadii P(ii,a) = 1 for every m, we 
conclude that Fo{a,u) = 1. 

Also, C5(u) appearing in F2, has to be zero. We proceed 
further in a similar way. In order to find Ci{ii) entering 
Fi{a,ii), we have to solve Eq. (PT|) for 71 = 3, and the 
procedure continues for higher-order terms in m. To be 
able to find Fi we have to solve the differential equations 
(|4T|) for al\n<i + 2. 

Apart from the already stated Eqs. (P5|) and (|44p we 
give the final expressions for 



F^ia, ii) = -ii^v In(-u) i—ct' + iaii + 603^^) 
6 

-lOa^ -t- 135a^w + IBOa^cs.ii^ 

^ 12960 

27a^M^(20u - 2Qv - 3) - IZbQh^c^v? 

12960 

45a3u3 (-4 (12c5 + h)ii + 2'iu'^ - 24^2 + SOw - 17) 



+ 



12960 
3240a2c3u'^(w-u) 

12960 
U5au* (-8 (6c5 -5)u + 2Au'^ ~ 2\xP- + 20i' + 9) 

12960 
270^5 (c3 (2MP- - 24i;2 + 36w - 5) - 48c7w) 



12960 



(A5) 



Higher orders can be calculated in the same manner. 

As an approximation for the tree model that seems 
justified for small enough mass, and large enough driv- 
ing velocity, one can write 1 = /„ da J_ dwP w 

/^ da /p°° dwP'^' . From this condition follow Eqs. ([9T|) . 
(|92]) . and we find 



r(i)) 2r(i>) 

X {-2iiv In (u) + 2mw?A(") [y) + u^ - w^ + v\ 



24r(D) 
X f — 12ww In (ii) [~iiv In (ii) + ii (ii + 2) — v^ + v) 

+ \2iiv\^(v) 

X (— 2Muln(M) -I- uv'ip{v) + ii{ii + 2) — v'^ + v) 

— iiv^^ '{v) ) 

- 6u'^{v - 4)w 4- 6u(5 - 3v)v + 3u'^ + lOu^ 

+ {v-2){v-l)v{3v + 5))] (A6) 

A similar analysis can be done for the yit model, and 
there Eq. ()A6P is obtained for small mass after neglecting 
complex velocities, as shown using the instanton solution 
in Sec. EB 
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Appendix B: Matching of it ~ m and u ^ m at small 

mass 

Next we discuss some remaining details of the pertur- 
bation theory presented in Sec. IIVBI Ahhough we could 
not solve analytically the equations that determine the 
solution in region 2, in this appendix we demonstrate that 
we properly organized the perturbation theory in region 
2. We will prove the matching between the distribution 
function in regions 1 and 2, without explicitly solving the 
equations in region 2. 

The matching condition at the boundary reads 
P(i)(w,a) sa P^'^\u,a). From that follows 



where c" = for i < as well as for i > [3n/2]. Indeed, 
our results satisfy this structure. 



Note that for i = we obtain Cq 3n 



?i-i 



L.0 72 = 0. To 
conclude, using Eq. (|B6p and knowing Cq = 1 wc can 
find all others coefficients cf. Then the distribution of 
velocities has for large enough velocities the form 



p^'\^ = T. 



■,v-k-l 



3k 



y2^'^-^+^/^ct''T{3k- 1 + 1/2). 

(B7) 

"(2) 

Now we consider the next term in the expansion P^ . 
Knowing P^^' . from the matching condition follows that 



n=0 

00 

+ In {m)m''-^ ^ ^,l^^(a, u)r 



n=0 



(Bl) 



where . . . denotes that there are other terms ^ (In to)* 
where i > 1, as well as ^ ■m"~^^'^ . Note that In (m) terms 
come from the ln(M) dependence in P^^\ while ~ m'"~^^'^ 

comes from terms like ^ C3 in P^^\ Note, that Vn satis- 

~ (2) 
fies the same equations as PA , but with different bound- 

ary conditions, i.e. it contains only some terms from P„ ' 

(see below). The same holds for the omitted terms which 

can be analyzed in the same manner. 

For large enough a, ii 



P^^\a,{L)^e 






Fn{a,u) ^ lim^^o 



^ T(v)V2ti 

n— v / ^ 

Fn {a^/m, iim) 



Fn{a,u), 



(B2) 



(B3) 



where F„ are determined by Eq. (jSJ. Analyzing expres- 
sions for Fn, we find 



F„{a,u) 



[3n/2] 
i=0 



(B4) 



where [x] rounds x to an integer such that [x] < x and c" 
are numbers. Strictly speaking Eq. (|B4p holds for n < 7 
and for n > 7 there might be some additional (Inu)' 
terms, but in that case the discussion would be similar. 
By plugging Eq. (JB1[) into Fokker-Planck Eq. ([M)) one 
finds 



dp;^ 



(2) 



d- 



)^r)-£.(-r)=o 



da 



In order that both Eqs. ft32|) and (|B5j) hold, 

cj'(3n - 2i) - <_i(3n - 2i + 2)(3n - 2i + 1) 



(B5) 



1 



:/(-:^ 



2n + z) + ^c^ 



+ c"r2^w(3n - 2i + 1) = , 



(B6) 



Pf (a,^)«e-t^g 



^T(v)V2n 



iia,u), 



(B8) 



[3n/2] 



Qnia,u) ^ — uFn{a,u) +u ^J «' 



3n-2^:r^^pn 



i=2 



[3r,/2] 



-t- u In u 2, 1^ 



Inu ^ a3"-2i^*d^, 

i=3 



(B9) 



for large enough velocities. On the other hand, the equa- 
tion for P{ (a,u) reads 



dPl 



(2) 



_9_ 

i(2: 



S(2) 



(« - a)Pi^ 



0. 



5u 

From Eqs. (|B8)) and (|B10|) one obtains 
p^S'^ - 2i) - P^-i{in -2i + 2)(3n - 2i + 1) 



(BIO) 



-^-(^ 



1 



2n + i) + ^r""' + flZ.^v{in - 2i - 

- c'/,-' (3n - 2i + 1) + d^7i' = 0, 
dn3n - 2i) - (i"_i(3n - 2i + 2)(3n - 2i + 1) 



1) 



(Bll) 



+ dU{ 



1 



2n + i) + -d" 



+ (flZ^'"{in-2i + l) = 0, 
(B12) 

where p" = for z < 2 or i > [3n/2], and d" = for i < 3 
or i > [3n/2]. Note that knowing d\ — v we can find all 
other coefficients. By examining the expressions for Qn 
that follow from P„, we found that both Eqs. (JBIII) and 
(jB12p hold. One can repeat the procedure for higher- 
order terms in the same way. 

Similarly, we find that in region 1 



,(2) 



0, 



n 



(2) 



r„_i_2„ [3"/2] 

T{v)V2^ S 



^'ii'd?. 



(B13) 
(B14) 



The analysis presented in this appendix confirms that 
perturbation theory in region 2 is properly organized and 
matching between the distributions in regions 1 and 2 
holds. 
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Appendix C: Large-deviation function 

Here we comment shortly on the Laplace transformed 
Fokker-Planck equation for the tree model and then 
present a calculation of its large deviation function. It 
allows to show that the equivalence of the large-deviation 
function for the tree and \/ii model for positive velocities 
holds self-consistently. It may allow to obtain a complete 
proof of the equivalence if one subcase (see below) could 
be ruled out, but we have not succeeded at this stage in 
doing so. 

The Laplace-transformed Fokker-Planck equation 
reads 



dP+ dP+ / K \ dP+ I K K^ 
dt Ok \ mJ d\ \ m m? 

K 



m 



-vP++<^iK), 



dP- dP- / K\ dP- f K K^ 



dt Ok V m/ dX \ m m" 

K 



m 



-vP- - <I>(k), 



(CI) 



(C2) 



where 



*('«) == ^ / dae'^"P(u = 0,a) (C3) 

OK J_oo 

Here we kept only the one boundary term at u = 0, 
assuming that other terms at "infinity" vanish. Here 
P+{X,k) = /g°°duX^dae^"+''''P(M,a). We define P_ 
similarly, with the difference that the integration is over 
negative velocities. 

Eq. (jCip (which describes the contribution coming 
from positive velocities) can be compared with the cor- 
responding equation for the v^ model ()108p . One sees 
that the main difference is the additional boundary term 
$(«;). When this boundary term vanishes, the two mod- 
els become equivalent. 

In the stationary case, introducing the time parame- 
terization as in Eqs. (jllip and ()112p we obtain 

^ - ^-h ~ ^i<t)) = 0. (C4) 

The solution of this equation is 
P+(Ao,.Jo)=i^+(0,0)e"^(^°'«°) 

(C5) 



where Z is given by Eq. ()118p with K{t) = kq and X{t) = 

Ao. 

Introducing F{x,y) as in Eq. (|194[) . we find $(k) = 
„^g«[G(ft)-/o]^ where P+{0,0) = e-'"f« = (.-vvaiu^,yF{x,y) ^ 
Here 

G{k) - /o = -minj,{P(0, y) - ny}, (C6) 



and ao is the corresponding value of the acceleration 
when the minimum is reached. If we assume that the 
first term in Eq. (|C5|) gives the main contribution in the 
limit w — >■ cxj, we obtain 

-/o + Z{X, k) = -min^^y{F{x, y) ~ ny - Ax}. (C7) 

Then indeed the second term is smaller than the first one 
in Eq. (|C5|) . since it reads 



-t-[Z(A(s),K(.s))-G(«(s))]^^ 



Q-vfopyZ{\,K) 



g-'"/Ogf^(A,K)g-'«min^£(_o„ t){Z(A(s),K(s))-G(K(s))} (Q,S) 

and using (jC6p and (|C7[) it follows that 0(Ao,ko) = 
min,g(_oo^t){Z(A(s),K(s)) - G{n{s))} > 0. 

However, if we assume that the second term gives the 
main contribution, then (jC8[) implies that 



- /o + Z{X,k) - 4>{X,k) 

= -mmx,y{F{x, y) - ny- Xx}. (C9) 

From this equation and Eq. (|C6p follows that (/)(Ao, kq) < 
Z(Ao, Ko) — G{i^o)- This statement is not in contradiction 
with the assumption, but on the other hand it does not 
follow from (|C9P that (/)(Ao,ko) < 0, as it should be if 
the second term is the dominant one. It would be nice 
to show that this possibility is ruled out, which would 
provide a proof of the equivalence of large deviation func- 
tions, independent of the one given in the main text. 

Note, that /o is expected to be zero and then (|C7p is 
equivalent to Eq. (|195|) . 



Appendix D: Moments of the distribution function 
for Vm model 



Using the equations derived in Sec. IV CI we find the 
moments characterizing the distribution function of the 
\Ai model discussed in Sec. IV CI Choosing t* = with 
a constraint u{t > 0) = and 11(0) = — A/ttt., we can 
rewrite Eq. (|115p as 



[mdt ~dt + l)u - w^ = XS{t) 

Then Z(0, A) = J°^ a{t)dt. 

One can write u = X^J^o """A", where 



(Dl) 



{mdf -dt + l)-u„(i) - ^ Ui{t)u„_e{t) = 6n.i6{t). 



e=o 



(D2) 



A solution of this set of equations is 

uq = 0, 

/n 
dTR{t - t) ^ ui{T)un-i{T) for n > 1, 

(D5) 



(D3) 

(D4) 
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with 

{md'i ^dt + l)Rt^t' = Sit - t'), (D6) 

^=e*/(2™) sini^V^^^T^), m > 1/4 



Rt 



V4^ 

^^gt/(2m) gjjjj^ (_1.^1_4j„)^ to < 1/4. 



(D7) 



The case m = 1/4 was treated in equation (jl47p . Using 

that P(A,k) = e''^^^'") and a5J5^P(A,K)|(o,o) = "^^o", 
we obtain the first few moments exactly: 



u = V, 

ii? = v{a + u), 
4^2 



m' + 2 
3 



+ 3crv + v' 



■A - 



U'^ = V 



f 60-3(5,71' + 6) cr2(3m' + 22)w ^ , 3\ 

H K + 2)(4m- + 3) + m^ + 2 + ^-^"^ + -> 

48cr4(TO'(103m' + 198) + 72) 10cr3(3im' + 30)v Str^ (3m' + 14)^2 



^ (m' + 2) (m' + 6) (4m' + 3) (9m' + 4) (m' + 2) (4m/ + 3) 
— _ ( 240cr5 [m' (m' (m' (69577i' + 4396) + 7666) + 4284) + 720) 
""v (m' + 2)2 (m' + 6) (4m' + 3) (9m' + 4) (16m,' + 5) 

2cr-*(m'(m'(9m'(225m' + 4138) + 130916) + 137592) + 39456)w 



lOcrw^ + v^ 



(m' + 2)2(m' + 6)(4m' + 3)(9?7i' + 4) 



15cr3(m'(4m' + 105) + 90)^2 5cr2(9m' + M)v^ 



+ 15crw'^ + -ty- 



(m' + 2)(4m' + 3) m' + 2 

^ /2880cr6(m'(m'(m'(m'(22015m' + 190244) + 521534) + 594996) + 277920) + 43200) 

II' - — ifX 



+ 



(m' + 2)2 (m' + 6)(m' + 12) (4m' + 3) (9m' + 4) (9m' + 10) (16m' + 5) 
168cr5(m'(m'(m' (20438m' + 114479) + 179748) + 93924) + 15120)d 

(m' + 2)2(m' + 6)(4m' + 3)(9m' + 4)(16m' + 5) 
14cr4 (to' (m' (3105m'2 + 39756m' + 125468) + 122256) + 33408) v^ 



(m' + 2)2 (m' + 6) (4m' + 3) (9m' + 4) 
105cr3(m'(4m' + 53) + 42)^3 35cr2(3m' + 10)w'' 



(m' + 2)(4m/ + 3) 



21crw^+w^ 



(D8) 
(D9) 

(DIO) 
(Dll) 
(D12) 



(D13) 



(D14) 



r 



where ra' — mii^/rf', a — ' / ^ .^J' . Here we used 
dimension-full quantities. These results are in agreement 
with the results from Section FVDI 



Appendix E: Perturbation theory for the instanton 
solution 



Here we give more details on the derivation of the per- 
turbative solution of Eq. (|115p considered in Sec. IV El We 
analyze the matching of the expansions in the different 
regions (|133p and (|137p at t ~ — m, in order to find ?/„ 
appearing in (|137l) . 



After determining /„ (see Sec. IV E[) we notice that for 
ri > it has the structure 



/„(a;) = 5]e'=-c(")(x,A) + y:4")(A)(-x)^ (El) 



fe=i 



k=0 



From Eq. (jl35p . and by comparing the coefficient in front 
of x^ we obtain 



{k + l)Al\'^ + 4^) - Y: AtP4' = 0. (E2) 

£=0 
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We find that 



1W/^^ _ ^k{^) 



Arw 



k\ 



5,(x) = ^a;"r 



This comes from 



n=0 



(E3) 
(E4) 



order terms j/„ can be found easily. We state here: 

6 
+ — e2*A (39A - I2t\ - BAX^ + 24tX'^) + 

-^e'Af - 72 + 36t - 6t^ + 192A - 144tA + 12<2A 

- 126A^ + 132^2 - 12t^A2 

+ ^A(72 + 36t + 6<2 - 231A - 114<A - ISt^A + 182A2 
+ 84^2 + 12^2 A^), (E8) 



k , . k / , \ oo oo 



=0 ^ ' £=0 ^ ' a=0 



b=0 



- ^ (a + 6)^t'^+'' = 5]t'x*(t + 1) = 5fc+i(.T) + 5fe(a;) 

a, 6=0 t=0 

(E5) 



By muhiplying Eq. (p5|) with 1/fc!, we obtain Eq. (|E2|) . 
if ^^^(A) = Sk{x)/k\. Using the "boundary condition" 
So{x) = 1/(1 - x) = yt|,°' = A, we find Eq. (|E3]) . 

It follows that the matching condition holds in zeroth 
order in m if 



2/0(2:) =: ^ S'n 



n=0 



A-l\ (-x)' 



A / n! 



A 



A + (1-A)e-^ 



(E6) 



2/2 (i) 



12(e-*(-l + A)-A)3 
X (\{ - 96 + 6t2(-l + A) + 141A - 25A2 

-6f(-8 + 9A)) 

+ e"*(72 + 6t2(-l + \f - 135A + 41A2 + 25A3 

-6i(4-9A + 5A2)) 

- 12(e"*(-l + A) + A)(-6 - 2f(-l + A) + 7A) 

X In [e^* + A - e"*A] + 24(-l + A)x 

X (e-*(-l + A) + A) In [e^* + A - e-'\f\ . (E9) 



Appendix F: The behavior of Z{\) for A -^ ±ioo in 
the x/w-model 

In this appendix, we derive the large-A behavior of 
Z{X) for the V^-model. 

In order to facilitate our thinking, we denote r := —t, 
such that the time r of the instanton goes from zero to 00. 
The behavior for A — > 00 is complicated. Analyzing Z{X) 
for real A, one finds that there is a branch-cut singular- 
ity. On the diagonal in the complex plane, A ~ 1 -I- i, the 
instanton solution of (jlOll) looks rather chaotic. For com- 
plex A with vanishing real part, the behavior is slightly 
simpler: In the complex plane, u(t), its derivative u(t), 
and the energy 



In order to have a smooth function u[t) at i ~ —m in 
each order ?i > in in, it should hold 



E{r):='^[dMr)f + ^u{rf--^ (Fl) 

are behaving as indicated on Fig. [22] (blue solid lines). 
The energy is dissipated as 



k=0 



(E7) 



Then, ?;„>o(0) == X"^ and we find Eq. (|TiT]) . Higher- 



-E{t) = ~lidrU{T)Y 



(F2) 



We are interested in the case when A is equal to the imag- 
inary unit times a large positive number. We start by ne- 
glecting dissipation and the term ^ /^^ . One can check 
later on the trajectories we find that this approximation 
is justified. 
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FIG. 22. (Color online) Evolution in the complex plane of «(r) (top left), u{t) (top right), i?kin(T") (bottom right), and 
J„ dsu(s)ds (bottom left), for A = 50 000i, m = 1/4. The movement of ?i(r) starts at 0, moving in counterclockwise circles 
inwards. The kinetic energy evolves from left to right, and J7 dsu{t)ds from zero to the left. The dashed lines are the solution 
for /i = r; = 0, given in (|F6|I and (JFTJI . 



We will be working in dimensionless units a = rj ~ ^ = 
1. We have to solve 

ma^u(i) =w(T)^ (F3) 

driiiO) = X/m, (F4) 

w(0) = 0. (F5) 

The solution for drU = it can be parameterized by 0, 



A 



u{(t>) = {2e">' - 1) — 



(F6) 



While this is of course always possible, we claim and 
check in Eq. (jFSp below, that cj) is real, and thus has the 
natural interpretation of an angle. Using energy conser- 
vation, E{{)) = ^, and E{t) = l^[dru{T)f - ^ , this 
yields 



"(0) 



-6A2 



-e 3 



</T 



(F7) 



This allows to obtain the "angular velocity" (attention: 
the prime indicates the derivative w.r.t. the argument 4>) 



M'(0) 



|Ae^^^^^(l-e^^)^/^ (F8) 



Note that for e [0, 27r], expression (|F8[) is real, thus the 
curves are indeed parameterized by a real angle 0. 
Noting the differential dt. 



dt 



this allows to calculate the period as 



(F9) 



T = 



22/335/6r(i)r(7^3 






5.56022 



A 



(FIO) 



Next we switch on dissipation. It leads to a change in 
energy, for a time T or the equivalent phase $, 



E{^) - E{0) = E{T) - E{0) = / E{T)dT 

il{TfdT = - ^(0)2 J^ , (Fll) 
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In principle, the exact trajectory has to be put. Using 
the dissipation-less one, we find for one period 

E{27:)-E{0) «-- 



-I)r(i] 



2v^v^0F^ 



A 



3.33613 



A 



m \iy/m 

5/3 



If we call n the number of periods, then 
A2 



EiO) = ~->0 



-i)r(i) 



dn^ ^ ^J 2^^^^/m 



2E{2iin) 



-3.33613 



2E{2Tm) 



5/6 



(F12) 
(F13) 

5/6 

(F14) 



where we have re-expressed A by the energy itself, in 
order to allow for an iteration. Defining Aofr by 

Aeff(2^n)2 



E{2Tm) 



2m 



(F15) 



we get 



d Aoff(27rn) d 
i dn 

r( 



dn 



^-2mE{2nn) 

- 2E{2TTn) 
Aoff(27rn)"'^'^ 



-I)r(i: 



2^^0F 



1/3 



rM)r(i) 

2^^y^ 



-3.33613 



Aeff(27rn) 



2/3 



(F16) 



The integral over one period of L dtu{t), which con- 
tributes to Z{X) is 

'■^dru(i)= r ^u{cb) 

Jo 0(0) 

1/3 



^r(-i)r(|) /A 



22/30F 



-3.91452 



Am 



771 

i 

1/3 



(F17) 



Therefore, using that Z{X) = L u{T)dT we have the 
following relations: 

1/3 



d?! 
This yields 



Z(A) = 



^r(-i)r(| ) fx 

22/30F 



(F18) 



d{X/i) 

2^V^r(-i 



dZ(A) 

dn 

dA/i 

dn 



) f A 



r(-i)r(i) 



-1.17337 



1777. 
1/3 



-1/3 



(F19) 



where the additional minus sign has been introduced due 
to the fact that we now integrate in the opposite direc- 
tion. Integrating over A, and using that Z{Q) ~ 0, we 
find the asymptotic behavior 



>(A) 



15\y6V^r(|) /Am 



r(^) 



2/3 



= -1.76006 



\ \ 2/3 

Am ^ ' 



(F20) 



As an example, for A = lO^z, m ~ 1/4, our formula 
give Zasymp(lO^i) = —150483.0 whereas numerics gives 
Z{10'^i) = -150483.6 -H 330.444i. The subleading imagi- 
nary part is consistent with 



^guess(A) = -1.76006 



/ Xm\ f Xra\ 

(-) ^^H-) 



1/3 



(F21) 

On the negative imaginary axis, the result is the same, 
i.e. on the whole imaginary axis 



-^asynip^'^ j 



i5^V^r(f] 



(|A|m)'/^ . (F22) 



Appendix G: Large m expansion to second order 

Here we state the second order of large m expansion 
of the distribution function for \/ii model, see Sec. IIVDI 



(Gl) 



Q2 = T^ir"" cos(60) + 48r\-2 sin(30) 
576u* 

-f 432™^ sin(6') (r^ - Av) 

-hl5r2 cos(26i) (r^ - 12w) (r^ - 4w) 

-f 10 {r^ - 72r2u2 -f- 96w^) + 6r^ cos(46') (r^ - lOw) ) . 

It can be rewritten as: 
1 



Q 



■(85^ - 75a\' + \%h^v{u - v){hu -f 3z;) 



144^4 

+3u(15u'' - 28m^u - Qi?v{v -h 20) + I2uv^{iv + 8) 
+w2((24-17i;)w-f-80))). (G2) 

Integrating out a, we find the distribution for the veloc- 
ity: 



P{u 



1 -i^-..)-^ 
-_e-^r. — (1 



1 



-{lbu^~2%u^v (G3) 
v27ri' ' 48mu3^ 

~Qi?v{v + 15) + I2uv^{iv + 7) + v^{{& - 17u)u + 45)) 
\ m^ 



Appendix H: Exit probability 

In this section we calculate the probability E{a, u, t) 
that a particle starting at i = with acceleration a and 
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velocity it > had a negative velocity at some moment 
before or at time t. We call it the exit probability. The 
calculation is valid for all the models studied in the pre- 
vious sections. 

The equation for the exit probability reads as: 



dE dE u d^E 1 , ^dE , ^ 

+ —i~u-a + v)— (HI) 



dt dii m? d^a m 



with the boundary conditions: 



E{a <0,0,i) = 1, 
E{a,u> 0,0) = 0. 



(H2) 
(H3) 



We assume that the particle starts with a finite accel- 
eration and velocity. Then, introducing E{K,X,t) = 
/_ da /„ duE{a, it, t) exp (na + Xii), we obtain 



dE _dE 
'dt~'dX 



K 

m 



dE 
dn 



-A + 



^(1 

m 



i + /(«)' 



(H4) 



for we used that k > 0. The last term 1/k'^ + /(k) comes 
from the boundary term at li = obtained by performing 
the Laplace transform with respect to ii. Here /(k) = 
— /p daaexp (Ka)£'(a,M = 0,i). 

Using the method of characteristics we find: 



A 



m^ m 

K 



k = \ 

m 

dE E{t) 



{l-VK{t)) 



1 



di m ^ ^ " i^'^i't) 

where we used K(t) > 0. Then, we obtain 

1 



(H5) 
(H6) 

/(«(f)) = 0, (H7) 



E{ko > 0,Ao,to > 0) = / ds 
Jo 



«2(,) 



fi<^)) 



1 



to 



X exp 1 — / (1 - VK{s'))ds' ] , 



where nit) satisfies the following equation 

,2 



kit) 



K 

m 



K 

in 







(H8) 



(H9) 



with conditions ^(to) ~ '^o a-nd ^(to) = Aq — kqIui. If we 
introduce ui^s) ~ —K.{—s)/m, then u{s) satisfies Eq. (jllSp 
where in Eqs. (|116lll7p it has to be made a change kq— > 
— kq and Aq— > — Ao. The equation should be solved self- 
consistently, since /(k) is determined by E{a > 0,0, i). 
Increasing the driving velocity v and decreasing the mass, 
the boundary term /(k) decreases. It is expected that 
the approximation /(k) = becomes reasonable good 
for sufficiently large driving velocity and small mass. 



Appendix I: Sketch of a proof for convergence of 
large-deviation function 

Suppose that we have a numerical simulation of one of 
the models discussed in this work, at driving velocity v, 
which gives N data points. This allows to estimate Fy (x) 
in a certain domain x- < x < x+, with F{x-) = F[x+). 

Let us first estimate the probability p that there are 
data points left of x^ : 

da;i;e-"-^"(^) w N- — ; ^ w jy^-vF^x-.) 

F'ix-) 

(11) 
This means that the simulation has to be done at veloc- 
ities large or equal to Vc in order to satisfy (|IT|) . with 



HN/p) 
F(x„) 



(12) 



where we have replaced Fy(x) by the limiting function 
F{x). (For N = 10^°, and p = 10"*, this would give 
Vc ~ 21.5; using for Fy{0) the function for the -v/u-niodcl.) 
Let us now estimate which domain of the function 
Fy{x) we can estimate with relative statistical error 
smaller than e, at this velocity Vc- The number n of 
events in the bin around x of size S is 



-"cF(x)^^ > 



1 



and must as written be larger than 1/e^ 
yields 



Fix) < 



ln{e^SN) 



Fi^o] 



(13) 
Solving for F{x) 

\n{e^5N) 

HN/p) 



Fix^) 



1 



InjeHp) 
In AT 



(14) 



where in the last line we have supposed that TV is large. 
This is probably a rather crude estimate, but shows that 
for TV — > c» one can estimate F{x) for all x for which 
F{x) < F{x^). 

Let us now consider the yu- model, for which Zy{X) 
does not depend on u, and for which Fy{x) converges 
against the large-deviation function F{x). The above 
shows that there exists a simulation, which can estimate, 
with any given precision, F(x) for all x with F{x) < F{0) 
(see plot [201) • For further reference set x_ = 0, and 
x+ = 3.47268 the other root for which F{x+) = F(0). 

We remind that with probability 1 — p w 1, the sim- 
ulation has never encountered a negative velocity. Now 
repeat the simulation with the same parameters, with one 
of the other models. With the same probability 1 — p, 
these simulations have no negative velocities, and since 
then the particle will only move forward, give the same 
trajectory, and thus the same large-deviation function, 
within the (small) error-bars estimated above. We have 
thus proven Eq. (|192p of the main text. 

The only circumstances where the above argument 
might go wrong, is if there are strong correlations in the 
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tails. If e.g. rare events are correlated such that when- 
ever one gets one rare event (of negative velocity) then 



one gets with higher probability another one (clustering 
of rare events). In this case even the existence of a large- 
deviation function may be questionable. 



[1 

[2; 
[3; 

[4] 
[5 

[6: 

[7] 
[8 

[9; 

[lo; 

[11 

[12 
[13 
[14 
[15 
[16 
[17 
[18 
[19 

[2o; 

[21 

[22; 
[23; 

[24 

[25 
[26; 

[27 

[28; 

[29; 

[30 

[31 



D. S. Fisher, Phys. Rep. 301, 113 (1998). 

S. Brazovskii and T. Nattermann, Adv. Phys. 53, 177 

(2004). 

P. Le Doussal and T. Gianiarchi, Phys. Rev. B 57, 11356 

(1998). 

F. Colaiori, Adv. Phys. 57, 287 (2008). 

S. Field, J. Witt, F. Nori, and X. Ling, Phys. Rev. Lett. 
74, 1206 (1995). 

G. Gruner, Density Waves in Solids (Addison- Wesley, 
Reading, MA, 1994.). 

O. Narayan and D. S. Fisher, Phys Rev. B 46, 11520 

(1992). 

F. Perruchot and et. al, Physica B 284-288, 1984 (2000). 

C. Reichhardt, Phys. Rev. Lett. 86, 4354 (2001). 

R. Chitra and T. Giamarchi, Eur. Phys. J. B 44, 455 

(2005). 

L. F. Cughandolo, T. Giamarchi, and P. Le Doussal, 

Phys. Rev. Lett. 96, 217203 (2006). 

D. Wilkinsion and J. F. Willemsen, J. Phys. A 16, 3365 
(1983). 

J. Stokes, D. Weitz, A. D. J. P. GoUub, M. Robbins, 

P. Chaikin, and H. Lindsay, J. Phys. A 57, 1718 (1986). 

S. Ramanathan and D. S. Fisher, Phys. Rev. B 58, 6026 

(1997). 

J. S. Urbach, R. C. Madison and J. T. Markert, 

|Phys. Rev. Lett. 75, 276 (1995) 

D.-H. Kim, S.-B. Choe, and S.-C. Shin, 



Phys. Rev. Lett. 90, 087203 (2003) 

B. Alessandro, C. Beatrice, G. Bertotti, and A. Mon- 

torsi, J. Appl. Phys. 68, 2901 (1990). 

B. Alessandro, C. Beatrice, G. Bertotti, and A. Mon- 

torsi, J. Appl. Phys. 68, 2908 (1990). 

P. Cizeau, S. Zapperi, G. Durin, and H. Stanley, 

[Phys. Rev. Lett. 79, 4669 (1997) , 

S! Zapperi, R Cizeau, G! Durin, and H. Stanley, 

[Phys. Rev. B 58, 6353 (1998} 

G. Durin and S. Zapperi, 

[Phys. Rev. Lett. 84, 4705 (2000)' 

P. Cizeau, S. Zapperi, G. Durin, and H. E. Stanley, Phys. 

Rev. Lett. 79, 4669 (1997). 

P. Le Doussal and K. Wiese, Phys. Rev. E 79, 051106 

(2009). 

P. Le Doussal and K. Wiese, Europhys. Lett. 97, 46004 

(2012). 

P. Le Doussal and K. Wiese, in preparation. 

A. Dobrinevski, P. Le Doussal, and K. J. Wiese, Phys. 

Rev. E 85, 031105 (2012). 

J. Schwarz and D. Fisher, Phys. Rev. Lett. 87, 096107 

(2001). 

J. M. Schwarz and D. S. Fisher, Phys. Rev. E 67, 021603 

(2003). 

A. Mehta, A. Mills, K. Dahmen, and J. Sethna, Phys. 

Rev. E 65, 046139 (2002). 

D. Spasojevic, S. Bukvic, S. Milosevic, and H. E. Stanley, 

Phys. Rev. E 54, 2531 (1996). 

G. Durin and S. Zapperi, J. Magn. Magn. Mater. 242, 



[32; 

[33 

[34 
[35| 

[36 

[37 

[38; 
[39 
[40 

[41 
[42; 
[43; 
[44 



[45' 
[46; 

[47 



[49; 

[50 
[51 

[52; 

[53; 

[54 



1085 (2002). 

G. Durin, F. Colaiori, C. Castellano, and S. Zapperi, J. 
Magn. Magn. Mater. 316, 436 (2007). 
S. Zapperi, C. Castellano, F. Colaiori, and G. Durin, 
Nature Physics 1, 46 (2007). 
W. Doering, Z. Naturforschung 3a, 373 (1948). 
V. Lecomte, S. E. Barnes, J.-P. Eckmann, and T. Gia- 
marchi, IPhys. Rev. B 80, 05441 3 (2009) 
D. Bernard and B. Doyon, arXiv , 1202.0239 (2012). 
L. S. Levitov, in Quantum Noise in Mesoscopic Physics, 
NATO Science Series, Vol. 97, edited by Y. V. Nazarov 
(Springer Netherlands, 2003) pp. 373-396. 

C. W. J. Beenakker and C. Schonenberger, Physics Today 
56, 37 (2003). 

P.-E. Roche, B. Derrida, and B. Dougot, Eur. Phys. J. 

B 43, 529 (2005). 

L Klich, in Quantum Noise in Mesoscopic Physics, 

NATO Science Series, edited by Y. V. Nazarov (Kluwer, 

Dordrecht, 2003). 

J. Gabelli and B. Reulet, 

■phys. Rev. B 80, 161203 (2009)', 

D. Gorokhov, D. Fisher, and G. Blatter, Phys. Rev. B 
66, 214203 (2002). 

T. Nattermann, T. Giamarchi, and P. Le Doussal, Phys. 
Rev. Lett. 91, 056603 (2003). 

V. Repain, M. Bauer, J.-P. Jamet, J. Ferre, 
A. Mougin, C. Chappert, and H. Bernas, 

Europhys. Lett. 68, 460 (2004) 
A. A. Middleton, Phys. Rev. Lett. 68, 670 (1992). 
The next order at small u should be examined with the 
same conclusion. 

In the limit v = 0'^ metastable states can still be defined. 
In the limit ri,v = 0^ one can define the quasi- 
static motion of the edges of the oscillation interval 
[umin(i), Umax(i)]. It should also proceed by jumps. 
J. M. Schwarz and R. Maimon, 

[Phys . Rev. E 64, 016120 (2001) 



[55] 



Y. G. Sinai, Theor. Math. Phys. 90, 219 (1992). 
T. W. Burkhardt, J. Phys. A: Math. Gen. 26, L1157 
(1993). 

P. Le Doussal and K. Wiese, Phys. Rev. E 79, 051105 
(2009). 

The function / does not enter, since the Fokker-Planck 
equation is a linear equation. 

Although this could be interpreted by saying that the ef- 
fect of the disorder in the moving frame is equivalent to 
an equilibrium white noise of variance 2'qTcff for a par- 
ticle in a quadratic well with damping 77 (i.e. satisfying 
usual fluctuation-dissipation relations), this would lead 
to an incorrect value for the noise; hence it is is not a 
valid interpretation. There is an underlying Hamiltonian 
system, but it is non-standard and involves the accelera- 
tion, see Sec. lVID] This is because the starting equation 
of motion can only be written in a Langevin form in terins 
of the velocity. 
D. L. Scharfetter and H. K. Gummel, IEEE Trans. Elec- 



42 



tron. Devices ED-16, 64 (1969). 

Note however that for m < there are solutions which 
vanish at infinity and do not satisfy (|102|) . 
The result is exact in the region of parameters discussed 
at the very beginning of this section. 
M. J. Ablowitz and A. Zeppetella, Bulletin of Mathemat- 
ical Biology 41, 835 (1979). 

A. D. Polyanin and V. F. Zaitsev, Handbook of exact 
solutions for ordinary differential equations (CRC Press, 
1995). 

F. Evers and A. D. Mirlin, Rev. Mod. Phys. 80, 1355 
(2008). 

Y. Fyodorov, J. Stat. Mech. , P07022 (2009). 
Note that here the natural unit of temperature is 



[56] 
[57] 
[58] 
[59] 

[60] 

[61] 

[62] 

ninvfi^ — n'^'/a'' as can be seen from (|217|) and using 
that the dimensionfuU instanton solution is UA'm'""(*) = 
iif^ -dimioss (t/r„). The model now depends on three di- 
mensionless parameters, m,v and T in the units defined 
here. 

[63] n here is noted m there. 



CONTENTS 



I. Introduction 



II. ABBM model 



Q 



III. 



ABBM model with inertia 

A. Definition of the model 

B. Phenomenology of the inertia model 

1. Qualitative features of trajectories 

2. Velocity distributions 



IV. 



V. 



Tree model 

A. Definition of the model 

B. Perturbation expansion in small m 

C. Large driving velocity u 3> f p 

D. Large mass m ^ 1 

1. D{u) = u model 

2. Tree model D{u) = |u| 

E. Numerical solution of the Fokker-Planck 
equation 

F. Comparison of the numerical solution of the 
Fokker-Planck equation with perturbation 
expansion in m and 1/m 

\/u model 

A. Strategy 

B. Definition and basic properties 

C. Connection between instanton approach and 
Fokker-Planck approach 

D. Moments of the distribution function 



M 



E. Perturbation at small m for P(u) 

F. Critical case m — 1/4 

1. Perturbative expansion for Z{X) 

2. Perturbative expansion for the critical 
instanton and Aj 

3. A good approximation for Z{\) from the 
branch cuts 

4. Large A-behavior 

5. From Z{X) to Papprox(u) 

6. The v^-model as an approximation for 
ABBM 

G. Exact results for m = 6/25 

1. Fisher-Kolmogorov equation 

2. Abel differential equation of the second kind 

VI. Large-deviation function 

A. Definitions and numerical determination 

B. Convergence of the large-deviation functions 

C. Large deviation function in perturbation theory 
for small m 

1. Expression obtained from previous results 

2. Equation from Fokker-Planck 

D. Large deviation function for large m and 
matching of small and large m 

E. Large-deviation function for m = 6/25 

VII. Thermal and quantum fluctuations 

A. Classical systems: Thermal fluctuations 

B. Quantum system 



VIII. Conclusions and discussion |32 

IX. Acknowledgments [33 

A. Perturbation theory at small m for P{u) for the tree 
model [33 

B. Matching of u ~ jn and u^ m at small mass [34 

C. Large-deviation function [35 

D. Moments of the distribution function for Vm model [35 

E. Perturbation theory for the instanton solution [36 

F. The behavior of Z(\) for A — > ztioo in the 
Vw-model [37 

G. Large m expansion to second order [39 

H. Exit probability [39 

I. Sketch of a proof for convergence of large-deviation 
function [40 

References [44 



